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ABSTRACT 


The  reliability  of  the  structure  can  be  calculated  by  combining  the 
non-uniform  stress  distribution  of  a  structure  from  finite  element  analysis 
with  the  statistical  strength  theory  under  two  dimensional  non- 
homogeneous,  uniaxial  stress  limitations.  The  specimen  models  with  and 
without  notch  are  the  sample  structures  to  illustrate  the  calculations.  The 
statistical  strength  of  the  structure  is  cast  in  the  standard  Weibull  form 

characterized  by  the  structural  scale  parameter,  $^  and  shape  parameter, 
<X|:  which  are  functions  of  the  material  scale  and  shape  parameters  (<x  and 

3.)  Therefore,  this  thesis  demonstrated  that  the  scale  parameter  of  the 
structure  remains  constant  for  different  load  magnitude,  only  the 
structural  geometry  and  loading  condition  cause  the  change  of  material 

scale  parameter,  3  to  the  structural  scale  parameter,   3^-  The  results  also 

shows  that  the  degree  of  uniformity  of  stress  within  structural  element 
and  the  linearity  of  the  load  magnitude  and  stress  distribution  within  the 
structure  affect  the  accuracy  for  the  calculation. 
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I.      INTRODUCTION 

A.     RELIABILITY  REQUIREMENT  FOR  COMPOSITE  STRUCTURES 

Reliability  theory  became  an  independent  scientific  discipline  at  the 
beginning  of  the  1950's  under  the  influence  of  the  rapid  development  of 
radio  electronics,  computer,  devices,  and  rocketry.  Modern  radio  electronic 
apparatus  and  digital  mechanics  consist  of  a  very  large  quantity  of 
components.  If  the  failure  of  one  component  leads  to  the  failure  of  the 
unit  as  whole,  then  evidently  the  possibility  of  high  reliability  (freedom 
from  failures)  will  diminish  rapidly  as  the  unit  becomes  more  complex.  In 
this  connection,  problems  arise  on  the  prediction  of  the  reliability  of 
planned  units,  on  the  development  of  measures  to  increase  the  reliability, 
on  a  fundamental  for  reliability  testing  methods,  etc.  All  these  questions 
are  the  subject  of  reliability  theory. 

The  mechanical  strength  and  stiffness  of  a  structure  are  one  of  the 
aspects  of  reliability.  The  engineering  design  concept  is  not  only  solving 
the  stresses  and  strains  which  originate  in  structures  under  various 
external  imposed  conditions  but  also  trying  to  make  them  to  be  operated 
sufficiently  reliably  throughout  their  established  utilization  time. 
Therefore,  in  the  concluding  stages  of  an  engineering  design  structural 
mechanics  inevitably  comes  into  resolution  with  reliability  theory.  For 
critical  application  of  composite  as  an  aircraft  structure,  the  study  of 
reliability  in  composite  structures  is  required. 
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B.      HISTORICAL  BACKGROUND 

Bullock  [  Ref.  1  ]  used  Weibull  theory  to  predict  the  strength  of  the 
composite  materials  between  different  laboratory  sample  configurations. 
He  investigated  the  strength  ratios  of  composite  materials  in  flexure  and  in 
tension.  He  assumed  the  shape  parameter  (the  exponent  characterize 
flaw-density  that  determines  the  scatter  of  strength  of  the  material)  was 
not  changed  when  the  geometry  or  size  of  the  material  was  changed.  He 
used  the  same  shape  parameter  for  the  tensile  test  and  flexure  test  which 
based  on  homogeneous  state  of  stress  and  heterogeneous  state  of  stress 
respectively.  He  found,  however,  based  on  the  experimental  results,  the 
hypothesis  that  the  shape  parameter  was  constant  was  not  realistic. 
Changes  in  shape  parameter  were  in  fact  observed.  Whitney  and  Knight,  n 
[  Ref.  2  ]  used  a  statistical  strength  theory  based  on  Weibull  distribution  to 
explain  the  difference  between  unidirectional  tensile  data  generated  from  a 
flexure  test  and  a  standard  tensile  coupon.  The  result  was  shown  a 
significantly  larger  variation  in  tensile  strength  versus  flexural  strength. 

Rosen  [  Ref.  3  1  presented  a  theoretical  and  experimental  treatment  of 
the  failure  of  composites,  consisting  of  a  matrix  stiffened  by  uniaxial ly 
oriented  fibers,  when  subjected  to  a  uniaxial  tensile  load  parallel  to  the 
fiber  direction.  He  observed  that  a  portion  of  the  fiber  at  each  end  was  not 
fully  effective  in  resisting  the  applied  load  because  of  the  axial  load  was 
transmitted  by  shear  through  the  matrix  to  adjacent  fibers.  Basically,  the 
fibers  failed  as  a  result  of  a  statistically  distributed  flaws  or 
imperfections  and  composites  failed  as  a  result  of  a  statistical 
accumulation  of  such  flaws  over  a  given  region.  Therefore,  the  ineffective 
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length  of  the  fiber  was  introduced  in  order  to  consider  the  composite  to  be 
composed  of  a  series  of  layers  of  this  length.  Then  the  segment  of  a  fiber 
within  a  layer  might  be  considered  as  a  link  in  the  chain  that  contribute  the 
fiber.  Each  layer  was  then  a  bundle  of  such  links,  and  the  composite  as 
whole  can  be  modeled  by  series  of  such  bundles.  Weibull  theory  was  used 
to  define  a  statistical  distribution  of  flaw  or  imperfections  that  result  in 
fiber  failure  under  applied  stress.  The  statistical  accumulation  of  such 
flaws  within  a  composite  material  was  demonstrated  to  be  cause  of 
composites  failure.  He  concluded  that  the  stress  concentrations  occurred 
in  fiber  adjacent  to  the  fibers  break  were  considered  only  to  balance  the 
effect  of  the  variation  of  ineffective  length  which  affect  to  the  composite 
strength.  Knight,  C.E.  [  Ref.  4]  used  a  finite  element  stress  analysis  and 
statistical  strength  theory  to  assess  the  influence  of  the  stress 
concentration  on  the  ultimate  strength  calculation  for  composites  from 
the  Split-D  Test.  He  observed  that  the  geometry  changed  causes  the 
nonlinearity  and  the  boundary  condition  would  be  changed. 

The  Probabilistic  Statistic  Failure  of  Composite  Materials  [  Ref.  5  ] 
demonstrated  how  to  develop  a  theory  that  combined  the  statistics  of 
composite  material  failure  with  the  orthotopic  nature  of  composite 
materials.  The  effects  of  loading  history  and  the  probabilistic  location  of 
the  failure  for  a  composite  can  be  accounted  using  this  theory.  In  order  to 
improve  the  efficiency  of  composite  designs  for  structural  analysis,  a 
theory  that  combining  the  statistics  of  the  experimental  failure  data  with 
an  orthotropic  material  failure  response  was  developed.  The  analytical 
theory    developed    extent    deterministic    failure    laws    to    include    the 
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statistical  scatter  observed  experimentally.  Tsai  and  Wu  [  Ref.  6  ]  are  one 
of  the  investigators  who  described  a  general  theory  for  the  failure  of 
composite  materials  as  the  strength  tensor  theory  but  determining  is 
complicated  by  wide  range  in  failure  properties  of  composite  materials  and 
scatter  in  experimental  results.  The  Weibull  weakest  link  theory  was 
developed  and  described  a  description  of  the  spatial  distribution  of  the 
failure  locations  with  a  sample  calculation  performed  for  the  case  of  a 
uniformly  loaded  tensile  specimen.  Three  point  bending  specimens  were 
used  to  predict  failure  location  by  applying  the  basic  weakest  link  theory. 
Four  point  bending  specimens  were  used  to  illustrate  the  effect  of  loading 
history  which  were  evaluated  by  applying  numerical  integration  to  the 
probability  of  failure  density  distribution  curves.  The  compressive  failure 
load  was  demonstrated  that  it  is  not  deterministic  failure  and  also  having 
some  scatter  which  is  independent  of  tensile  scatter  so  that  the  Weibull 
failure  distribution  for  tension  and  compression  are  different. 

In  ordinary  structures  having  many  stress  components  active,  a 
decision  on  the  probability  for  survival  must  be  based  on  statistical 
descriptions  that  admit  multiaxial  stress.  The  development  of  combining 
weakest  link  theory  with  the  strength  tensor  representation  is  required. 

C.      OBJECTIVE  AND  SCOPE 

The  driving  motivation  for  this  research  is  that  many  critical 
components  of  modern  aircraft  structure  are  made  of  composite  materials. 
However,  the  reliability  design  is  necessary  for  engineering  designs 
requires  not  only  the  stress  analysis  but  also  statistical  theory  of  failure 
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strength  of  composites,  The  latter  motivated  the  investigation  of 
statistical  theory  for  composites  failure  strength  where  are  summarized  in 
the  last  section.  Any  aircraft  structural  components  inevitability  have  the 
non-uniform  stress  in  the  neighborhood  of  edge-notches  and  holes.  Size 
effect  is  the  one  that  is  importance  for  the  composite  study  because  of  the 
difficulty  and  the  economy  that  one  could  not  experiment  the  several  large 
samples  to  predict  the  reliability  of  the  components.  Even  the 
experimental  test  samples  can  be  made,  the  difficulty  to  identify  the 
failure  location  is  also  impossible.  Finite  element  method  and  the 
statistical  theory  can  be  combined  to  explore  the  possibility  of  prediction 
both  the  reliability  of  the  structure  and  the  failure  location.  Basically, 
finite  element  method  can  be  used  to  calculate  the  stress  distribution  that 
occurs  in  the  structure  under  an  arbitrary  loading  boundary  condition.  The 
structure  is  divided  into  elements  be  sufficiently  small  that  the  stress  is 
uniform.  This  solution  of  physically  continuous  system  can  then  be 
replaced  by  the  solution  of  the  discrete  system  of  the  elements.  The 
statistic  Weibull  theory  can  be  used  to  model  composite  materials 
composing  of  a  series  of  bundles  which  are  in  turn  composited  series  of 
fiber  links.  Another  assumption  for  the  statistical  Weibull  theory  is  the 
local  stress  in  each  link  must  be  uniform.  For  the  connection  to  the  two 
dimensional  finite  element  method,  instead  of  using  link,  the  element  will 
be  introduced.  For  one  dimensional  case,  the  link  can  be  used  but  in  two 
dimensional  case,  the  element  has  a  width.  The  operation  of  size  effect 
formulation  based  on  Weibull  theory  and  finite  element  method  will  be 
presented  in  next  chapter. 
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The  scope  of  this  research  limits  to  the  treatment  of  statistical 
strength  of  composite  structures  which  are  two  dimensional,  subjected  to 
nonhomogeneous  uniaxial  stress  including  size  effect.  The  specimen 
models  are  plate  with  and  without  notch  subjected  to  the  tensile 
displacement  loading  case.  The  finite  element  code  name  NIKE2D  [  Ref.  7  ] 
is  used.  N1KE2D  is  a  two  dimensional  finite  element  program  which  is 
developed  by  Hallquist,  J.O.  [  Ref.  7  ].  It  can  be  used  to  calculate  linear  and 
non-linear  finite  element  problems.  In  this  thesis,  only  the  linear  portion 
is  utilized.  MAZE  [  REf.  8  1  computer  program  is  used  to  generate  specimen 
model,  grid,  define  boundary  conditions,  and  also  create  the  input  file  that 
can  be  immediately  use  for  NIKE2D  program.  ORION  [Ref.  9]  is  the  post 
processing  of  N1KE2D  which  can  be  used  to  create  contour  curves  of  stress, 
displacement,  etc.  The  observation  of  uniform  stress  in  each  element  will 
be  made.  The  observation  of  using  half  of  the  entire  model  by  symmetrical 
will  be  made,  which  bases  on  the  idea  that  the  model  can  be  divided  into 
smaller  elements  when  compare  to  the  entire  model.  The  linear 
relationship  between  external  loads  and  the  local  stress  in  each  element 
will  be  verified.  The  plate  without  notch  is  used  as  bench  mark  to  compare 
the  reliability  results  to  the  plate  with  notch.  The  specimen  model 
dimension  bases  on  specimen  dimension  which  using  in  the  shear  test 
experiment.  The  shear  model  that  involves  with  the  multiaxial  stress 
which  can  not  yet  be  treated  in  this  thesis.  Future  work  requires 
generalization  of  current  results  to  multiaxial  stress.  When  such  extension 
is  made,  realistic  prediction  of  composite  structures  will  then  be  possible. 
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II.      STATISTICAL  MODEL  OPERATION 

A.     REVIEW  OF  WEIBULL  WEAKEST  LINK  FAILURE  THEORY 

consider  an  infinitesimal  volume  element  SV  at  point  x  and  subjected 

to  stress  state   Ojj.   The  probability  of  failure   Sf   of  volume  element   SV 

is  assumed  to  be  given  by: 

5f  =  *  [  ay  [x]  ]  8V  (  2.1  ) 

where:  *  =  failure  parameter  which  is  a  function  of  Oy.  which  is  a 

function  of  position,  x. 

It  is  assumed  there  are  only  two  possible  states  for  the  volume 
element,  SV  ,  either  the  volume  element  survives  the  loads  or  fails.  From 
this  assumption  the  probability  for  survival  or  reliability  is  simply: 

SR  =  1  -  Sf  =  1  -  *[  Ojj[x]]SV  (  2.2  ) 

The  basic  assumption  in  weakest  link  theory  is:  If  any  volume,  SV  ,  in 
the  total  structure  volume,  V  ,  fails,  then  the  entire  structure  occupying 
volume  V  fails. 

Consider  the  structural  volume,     V,  to  be  divided  into     N     volume 

elements    SVk  ,  with  each  of  the  volume  elements  located  at  point    xk  . 
Then  the  total  volume  is  given  by: 

V  =    lim  Z  SVk  =  Jv  dV;         k  =  I,  2  .  .  .  N  (  2.3  ) 

where:  N  -  <*> 
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The  probability  for  the  volume   V    to  survive,   R  ,  under  weakest  link 
theory,  is  given  by: 

R  =  lim  TT  SRk  =  lim  TT  { 1  -  *[  Ojj  [xkl  ]  8Vk  }  (  2.4  ) 

where:  k  =  1,  2  . . .  N ,  and 

Equation   2.4  can  be  equivalently  written  as: 

lnR=   lim  Z   In  {  !  -  ¥[  Oj;[xk]SVk  }  (2.5) 

where:  k  =  1,  2  . . .  N  ,  and 

N-oo 

The  quantity    *.&Vk  is  small  and  from  the  Taylor  series  expansion 

In  (1  +  x)  ^  x  for  x  «  1 

Equation  2.5  becomes: 

In  R  =  lim  I  -{  *[o-j;[xk]]SVk}  (2.6) 

or  R  =  exp  -{  Jv  *[  Oy  [xk]  ]  dV  }  (  2.7  ) 

where:  k  =  1,  2  . . .  N  ,  and 

N  "►<» 

Equation    2.7  is  the  basic  equation  derived  under  the  assumption  of 

weakest  link  theory.  Note  that  the  stress  distribution  need  not  to  be 
uniform;  the  material  can  be  generally  anisotropic  and  the  material  can  be 
non-homogeneous.  The  equation  is  limited  in  that  it  describes  only  the 
reliability  of  the  entire  structure  and  provides  only  a  quantitative 
description  of  where  the  structure  will  fail.   The  reliability  becomes: 


R=  exp-fOjj/p)" 
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where:  o<  =  shape  parameter,  and 

0  =  scale  parameter 
In   the   next   section,   basic   theory   will   be   used   to   obtain   the 
formulation  of   size-effect   for   one   dimensional,      non-homogeneous, 
uniaxial   stress. 

B.      SIZE-EFFECT    FOR    NON-HOMOGENEOUS    STRESS,    ONE    DIMENSIONAL, 
UNIAXIAL  STRESS 

The  value  of     o<Lj   and     3Li      are  estimated  by  the  experimental 
procedure  from  the  given   n   numbers  of  sample  length    LI .   For  a  sample 

length   L  ,  the  evaluation  of    <xL  and   £L  can  be  implemented  according  to 

the  following  procedure. 

There  are     j     links  in  length     L  (Figure  2.1-a)  each  links  are  not 
necessarily  of  equal  length,  but  the  stress  in  each  link  is  homogeneous, ie 

the  stress  axp  within  the  pth  link  satisfies  the  following  condition: 

xoxP-{/lp  xaxdx}/ip  $  6  (2.8) 

where:  e  can  be  any  arbitrarily  small  value 

L  =  length  of  pth  link 

P  =  I,  2  .  . .  j 
The  reliability  for  each  link  is: 

Rp(o-xP)  =  I  -  f(axP)  =  exp  -(axP/flp)«p  (  2.9  ) 
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Figure  2.1       (a)     One  Dimensional  Links  Configuration 

(b)     Two  Dimensional  Specimen  Configuration 
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There  are  (lp/l)   unit  links  within   lp  ,  the  reliability  of  each  unit  link 

within    L  is: 

Rp(c/)  =  exp  -(<*XP/*  )"  (2.10   ) 

where:  ex.  and  3  for  unit  link  are  material  constants 

The  reliability  for  each  link   L  is: 

Rp(oxp)   =  TT  {exp  -(axP/3)°<}m  ;  m  =  1,  2  .  .  .    (lp/l) 

=  exp  {-(lp/l)(d  XP/P)°0 

where:  (L/1)  is  the  number  of  unit  link  in   lp 

The  reliability  for  then  entire  length   L  is: 
RL  =  R,  (a, )  R2(o2)  R3(a3)  R4(a4) . . .  Rj(Oj) 

=  exp  {-(li/1)(tf  i/W0")  exp  {-(l2/1)(o  2/3)oc).  •  •  exp  {-(Ij/1)(d  j/3)oc] 
=  exp  -{(lt/l)(o  i/3)(X  ♦  (I2/1X0  2/3)<x  ♦ . . .  ♦  (lj/0(a  j/3)<x) 
RL  =  exp  -{  2  l(lp/l)(o  xP/3n  )  ;  p  =  1,  2  . . .  j  (  2.11  ) 

Under  uniform  stress  condition  within  an  element,  let     op  is  the 

element  stress  which  is  constant  within    lp      length  element  and    dp  is 

calculable  from  a  given  boundary  condition  (B.C.)  characterizable  by  {P}, 
(eg,  a  load  vector).  The  stress  due  to  the  same  B.C.  with  a  different 
magnitude   P  is: 

op[P]  =  op[{P}]  (P/{P})  (  2.12  ) 

where:  p  can  be  of  any  arbitrary  magnitude  in  linear  system,  or 

P   is  of  neighboring  magnitude  of  local  linear  system. 
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tfp[{P}]  is  expressible  from  structural  mechanics  calculations, 
(eg,  a  finite  element  method.) 

In  order  to  reduce  the  value  of    oxp  which  is  a  function  of    P  in 

equation    2.11  to  a  single  random  variable  {P}  ,  the  substitution  of    tfxp[Pl 

by  oxI)[{P}]    is  made  within  the  same  manner  as  the  expression  in  equation 

2.12.  Thus,  the  equation    2.11    becomes: 

RL  =  exp  -{  E  [(lp/lXa  ^{PMPJHP/W*] }  ;  p  =  1,  2  . .  .  j 

Expressing  the  above  equation  into  the  Standard  Weibull  form: 
RL  =  exp  -{ (P/0L)°1 } 

Then:     06  =  ex 

3L  =  0/  {  2  [(lp/l)(a  xP[{P}]/{P})a]  },/a      ;  p  =  1,  2  .  .  .  j 

Let:  F  =  {  I  [(lp/l)  !/« (0XP[{P}]/{P})P  }-,/o<   ;  p  =  1,  2  .  .  .  j 

Thus:      3L  =  3F 

That  is,  given  a  structure  of  length  L  subjected  to  boundary  condition 
parameter  {P}  (  which  may  be  of  a  single  or  multiple  load  Pi  ,  P2,  .  .  . 
vectors  characterizable  by  a  single  parameter  {P}  )  which  give  rise  to  a 
unidirectional  stress   Op  computable  from  the  methods  of  mechanics  (  eg, 

finite  element  analysis  ).  L  are  element  lengths  which  are  segmented  as 
small  as  necessary  to  assume  the  stress  to  be  uniform.  Upon  sampling  n 
structures  measuring   Pn  failure  loads. 
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<xj_ ,  £[_  are  shape  and  scale  parameters  ( in  dimension  of  load 

estimated  from  D(Pn)  and  structural  size  L  . 

o<  ,  3  are  shape  and  scale  parameters  ( in  dimension  of  stress  which 
are  geometric  independent  material  parameters  ). 

F  is  a  conversion  factor  calculable  for  structural  (  finite  element  ) 
analysis. 

For  special  case,  homogeneous,  one  dimensional  uniaxial  stress: 

From:    F  =  {Z[(lp/l)1/o<(axP[{P}]/{P})]0<r,/o<  ;  p  =  1.  2  . .  .  j 
For  oxp  =  a  and  {P}  =  Pn  : 

F  =  [  (o/Pn)<*  {  Z  (lp/l)  ,/o<  }<*  r1/o<  ;  p  =  I,  2  .  .  .  j 

=  [  (o7Pn)  {  I  (lp/1) ,/a  }  f1  ;p=l,2...j 

=  (o7Pn)-' (l/l)"17* 

Thus: 

RL  =  exp  -{ [  (P/W/Co/P^"1  /(L/l)  ",/o<  ]<* ) 

=  exp  -{ (L/l)  [  (a/fl)(P/Pn)  r ) 

This  recovers  the  well  known  classical  result  demonstrates  the 
consistency  of  the  current  formulation  is  true. 

For  each  given  non-homogeneous  boundary  condition,  a  corresponding 
correction  factor  F  can  be  calculated  by  finite  element  such  that  the  size 
effect  can  be  predicted.  By  weakest  link  in  length  X  ;  it  is  inferred  that 
the  plane  of  facture  is  perpendicular  to  X  .  That  means  the  failure  plane 
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has  no  thickness.  Since  the  one  dimension  specimen  has  no  width,  fracture 
width  is  not  a  parameter. 

C.      SIZE-EFFECT    FOR    NON-HOMOGENEOUS    STRESS,    TWO    DIMENSIONAL, 
UNIAXIAL  STRESS 

This  is  the  extension  of  one  dimensional,  uniaxial   stress  to  two 

dimensional  uniaxial  stress.     The  definition  that  uniaxial  stress  implies 

failure  has  zero  dimension  in   X  direction  is  adopted  for  one  dimensional, 

uniaxial  stress  ,  which  has  been  described  in  the  previous  section.  For  two 

dimensional,  uniaxial  stress,  failure  automatically  requires  definition  of 

the  failure  plane  because  the  specimen  has  dimension  in     Y  direction 

(Figure     2.1-b),  and  it  is  no  longer  conceptually  a  point  as  in  the  one 

dimensional  case.    In  the  similarly  procedure  as  in  one  dimensional  case, 

the  specimen  is  divided  into  sufficiently  small  areas  within  which  the 

stress  in  direction  of  force  is  assumed  to  be  uniform.    It  is  seen  that  not 

only  the  length  of  the  specimen  is  divided  into  small  element  lengths  but 

also  the  width  must  be  divided  into  small  element  widths  in  order  to 

provide  element  areas  which  are  satisfies  uniform  stress  assumption.  Thus: 

xaxP-|pJxaxdx/lp*€x  foran   'p 

where:  p  =  1,  2  . . .   j 

y0xP  "  wqJ  yaxdy/wq  *  ey  foral1   wq 

where:  q  =1,  2  .  . .  k 

Weakest  link  in  X  direction  is  defined  by  Eq.  2.8  and  unit  length  in  X 
direction  is  defined  by  Eq.  2.9. 
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Weakest  link  in  Y  direction  is: 

Rq[0xp<*l  =  exp  -foxW/3pq}«pq 

where:  oxp(J  is  stress  in  the  qth  element,  and  q  =  1,  2  . . .  k 

This  means  that  if  any  one  of  the  qth  fail,  all  elements  fail.  However, 
the  unit  (metric)  width  can  not  be  arbitrarily  defined.  It  is  visualized  that 
within  each  metric  area,  there  is  one  domineering  crack  (largest  crack 

perpendicular  to    ox).    When  subjected  to    Oxp(^,  the  domineering  crack 
extends  without  bound.    From  fracture  mechanics: 

where:  kc  is  a  material  constant  and  can  be  independently  measured. 

ac  is  weakest  link  in  Y  direction,  which  can  be  used  as  the  unit 
metric  width. 

Thus:     acW  =  (k^a/^)2  for  element  pqth 

As  well  as  the  division  of  length  element  L  into  unit  length,  the 

division  of  width  element  is  necessary.  In  stead  of  unit  length   I  ,   acpcl  is 

used  in  order  to  indicate  the  critical  width  in  each  element,     thus,  the 
reliability  of  entire  element  is: 

Rpq[a„P<l]  =  exp  -{ (lpq/l)(w  pq/acPq)(axP<V3)°< } 

=  exp  -{ (lpq/l)[w  pq/(kc/cxPq)2](oxP<V3)0< } 
=  exp  -{ (lpq/l)(w  pq/kc2)(axPq)2*o</3« } 
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where:  pq  =  I,  2  .  .  .  jk  (  subscript  is  changed  from  two  dimensional 

array  into  one  dimensional  array.) 

The  reliability  of  entire  specimen  (which  is  the  structure)  is: 

RE  =  R,  R2  R3  . . .  Rjk 

=  exp  -{  I  (lpq/l)(w  pq/K2)(axPq)2+°730<  } .  pq  =  1,  2  . . .  jk 
From  the  linear  relationship  between  external  load  and  local  stresses: 
a„W[P]  =  axP<J[{P}]  (P/{P)) 

where:  o*xp^[{P}]  is  element  stress  which  is  calculated  under  certain 

loading  boundary  condition  while  the  magnitude  of  load  is  {P}. 

axpcl[P]    is    element    stress    which   subjected    to    the    same 

boundary  condition  as  above  but  different  in  magnitude,  the  magnitude  is  P. 

Thus:     RE  =  exp  -{  £  (lpq/0(w  Qq/kc2)((P/{P))o^[{P)])2+0</^  ) 

where:  pq  =  1,  2  .  .  .  jk 

Express  it  in  term  of  the  Standard  Weibull  for  the  structure,  thus: 

RE  =  exp -{ (P%)*E  (  2.13  ) 

where:    o<E  =  <x  ♦  2 

3E  =  {P}0(o</o<+2)  {  I  (lpq/l)(w  pq/kc2)(ox^[{P}])°<+2  }HA«*2)) 

where:  pq  =  I,  2  .  .  .  jk  (  2.14  ) 

The  values  of  o<  and  3  of  material  are  altered  by  the  boundary  condition 
and  the  geometry  of  the  specimen  (or  structure)  changed  but  they  do  not 
depend  on  the  changed  of  the  magnitude  of  the  load.  Therefore,  the  shape 
parameter,  <x  and  the  scale  parameter,  3  can  be  solved  for  any  given 
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structural  geometry  and  boundary  condition  The  finite  element  method  can 
be  used  to  calculate  the  stress  distribution  of  the  specimen  by  dividing  the 
specimen  into  small  elements  within  which  the  stress  element  is  uniform. 

For  the  simple  limiting  2-D  case,  the  specimen  without  notch,  the 
stress  distribution  for  entire  specimen  is  uniform.   Thus: 

oxp<J  =  o*x  =  constant 

0E  =  ({P}/ax)  0(«/«*2^c(2/«*2)  {  z  [lpqW^]  j(-l/oc*2) 

=  {P}  kc(2/o<E)  s(«/otVaxA(1/(XE) 
where:  A   =  E  [  lpqwpq  ]  ,  pq  =  1,  2  . . .  jk 

<xp  =  <x+2 
Thus:       RE  =  exp  -{ (A/kc2)  (P/{P})  (ax/ J0*)  J^E  (  2.15  ) 
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III.      APPROACH  TO  THE  PROBLEM 

A.  SPECIMEN  MODEL  CONFIGURATION  AND  COORDINATE  SYSTEM 

The  specimen  model  has  L99  "  in  length  and  0.968"  in  width  is  align 
in  the  y-z  plane,  (cartesian  right  hand  coordinate  system  is  used).  The 
configuration  and  the  coordinate  of  the  specimen  without  notched  and  the 
plate  with  90  degree  notch  configuration  are  shown  in  Figure  3.1.  The 
fiber  direction  which    makes  the  orthotopic  material  having  large  elastic 

modulus  (  E^  )  is  considered  to  be  along  the  y-axis  (Figure  3.2). 

B.  FINITE  ELEMENT  PROCEDURE 
1.      Material  Constants 

The  finite  element  program  requires  defining  the  material 
constants  [  Ref.  7  ]  for  orthotopic  elastic  to  be  used  in  stiffness  matrix 
calculation.  The  material  constants  of  typical  graphite  epoxy  composite 
which  used  in  this  thesis  are: 

Ejj    =    18.65    Msi    (Mega  Pound  per  square  inch) 

E22  =  «-21     Msi 

E33=  1.21     Msi 

u2]  =  0.0339 

u3|   =  0.0339 


u32  =       0.4 
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Figure  3.2       Specimen  Models  Coordinate  System 
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G12  =  0.84  Msi 

2.  Boundary  Conditions 

The  specimen  model  is  subjected  to  the  displacement  boundary 
condition  (B.C.)  in  order  to  be  able  to  use  the  entire  specimen  model  in  the 
calculation.    An  original  problem  the  specimen  model  is  subjected  to  the 

displacement  B.C.  value     Uj    (Figure    3.3)  on  the  right  hand  side  in  the 

positive  y-direction  and  the  left  hand  side  is  subjected  to  the  displacement 

B.C.  value    Uj   in  opposite  direction.    In  order  to  fix  the  specimen  model 

from  rigid  body  displacement,  the  left  side  of  the  model  must  be  constraint 
in  the  direction  of  load.    Therefore,  the  displacement  value  on  the  right 

side  of  the  model  becomes  two  time  of     Uj .     The  transform  of  the 

displacement  B.C.  is  shown  in  Figure   3.4. 

3.  Data  Planning  and  Purposes 

N1KE2D  is  two  dimensional  axis  symmetry  finite  element  code 
that  the  thickness  is  assumed  to  be  infinite.  The  plane  strain  analysis  is 
used  by  assuming  no  displacement  along  x-axis  and  no  rotation  around  y  and 
z  axis,  which  it  is  satisfy  for  this  purpose.  There  are  eight  finite  element 
program  data  output  to  be  discussed  in  this  thesis: 

3.1  Specimen  model  without  notch  subjected  to  the 
displacement  B.C.  value  0.01" 

3.2  Same  specimen  as  3.1  but  the  displacement  B.C.  value  is 
0.001".  The  purpose  of  these  two  cases  is  to  demonstrate  that  there  is  no 
stress  concentration  along  displacement  direction  and  the  linear 
relationship  between  displacement  and  local  stress. 
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Figure  3.3      Displacement  Boundary  Condition 
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3.3  Specimen  model  with  notch  subjected  to  the  displacement 
B.C.  value  0.01"  and  the  stress  calculation  bases  on  four-point  integration 
(Ref.  7  )  in  order  to  observe  the  uniform  stress  in  each  element. 

3.4  The  same  as  3.3  but  the  calculation  bases  on  one-point 
integration  at  the  centroid  of  each  element. 

3.5  Half  of  the  specimen  model  with  notch  subjected  to  the 
displacement  B.C.  value  0.01"  and  four-point  integration  is  used.  This  case 
will  be  used  to  compare  with  the  entire  specimen  model  in  order  to  insure 
that  the  half  specimen  can  be  used  because  of  the  symmetrical 
configuration. 

3.6  Same  as  3.5,  but  one-point  integration  is  used. 

3.7  Same  as  3.6,  but  the  displacement  B.C.  value  is  0.004. 

3.8  Same  as  3.6,  but  the  displacement  B.C.  value  is  0.001. 
One-point  integration  results  will  be  used  to  be  input  data  for 

statistical  program. 

C.       STATISTICAL  PROCEDURE 

After  obtaining  the  results  from  the  finite  element  program, 
four-point  integration  results  have  to  be  considered  to  insure  that  the 
tensile  stress  in  each  element  is  uniform  or  almost  uniform.  If  they  are 
not  uniform,  the  size  of  the  element  need  to  be  reduced.  However,  due  to 
computer  memory  limitation,  absolute  uniform  stress  may  not  be 
attainable.  One  should  realize  that  what  spatial  domain  that  the  element 
size  should  be  reduced  and  what  area  the  element  need  not  be  changed.  The 
symmetrical  of  the  model  can  be  considered  to  be  used  only  half  of  the 
entire  model  in  order  to  increase  the  number  of  elements. 
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The  statistical  size  effect  in  two  dimensional,  non-homogeneous, 
uniaxial  stress  program  in  fortran  77  is  written  to  be  used  to  calculate  the 
reliability  of  the  specimen  models  by  using  the  finite  element  results  and 
shape  parameter,  «  and  scale  parameter,  3  from  several  numbers  of  certain 
sample  size.  There  is  one  program  including:  I)  the  specimen  model 
without  notch,  2)  the  entire  specimen  model  with  notch,  and  3)  the  half 
(top  side)  specimen  model  with  notch. 

This  program  is  presented  in  Appendix  A.  There  are  12  output  data 
files  which  are  shown  in  Table  4.3. 
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IV.     DISCUSSION  AND  CONCLUSION 

A.     EXAMINATION  FOR  UNIFORMITY  OF  UNIFORM  STRESS  WITHIN 
ELEMENT 

The  four-point  integration  are  obtained  in  entire  and  half  specimen 

models  subjected  to  the  displacement  boundary  condition  (B.C.)  value  0.01. 

For  the  full  specimen  model,  the  element  number  544  (shade  in  Figure  4.1) 

is  considered,  the  tensile  stress  values  can  be  read  in  figure  4.2  .where  the 

stresses  of  this  element  are  tabulated  under  rows  2173,  2174,  2175,  and 

2176  (relatable  to  the  element  by  dividing  the  last  number  by  4  recovering 

544)  and  the  tensile  stresses  are  0.03663  Msi,  0.03738  Msi,  0.01933  Msi, 

and  0.01851  Msi  (Figure     4.2)  respectively,  which  demonstrated  that  the 

stress  within  element  is  not  exactly  uniform  but  the  differences  are  small. 

The  element  number  of  the  half  specimen  model  at  the  same  position  as 

element  number  544  of  the  full  specimen  is  442   (shade  in  Figure   4.3)  and 

the  tensile  stress  can  be  read  in  Figure    4.4,  where  the  row  numbers  are 

1765,    1766,   1767,   and   1768   and  the   tensile  stresses   are   0.03667   Msi, 

0.03742  Msi,  0.01937  Msi,  and  0.01861  Msi  respectively,  again  they  are  not 

exactly  uniform  but  the  differences  are  small.    The  element  number  192  in 

Figure    4.5  (full  specimen  model)  is  outside  the  notch  region,  the  tensile 

stress  values  can  be  read  in  Figure    4.6,  where  the  row  numbers  are  765, 

766,  767,  and  768  and  the  tensile  stresses  are  0.06891  Msi,  0.06888  Msi, 

0.06852  Msi,  and  0.06855  Msi  respectively  ,  which  shown  that  the  stress 

within  element  is  not  uniform  but  the  differences  are  more  uniform  than 

element  in  the  notch  region.     The  element  number  of  the  half  specimen 
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model  at  the  same  position  as  element  number  192  of  the  full  specimen  is 
230  (Figure  4.7)  and  the  tensile  stress  can  be  read  in  Figure  4.8,  where 
the  row  numbers  are  917,  918,  919,  and  920  and  the  tensile  stresses  are 
0.06906  Msi,  0.06905  flsi,  0.06871  Msi,  and  0.06872  Msi,  while  they  are  not 
exactly  uniform  but  the  differences  are  also  smaller  than  the  element  in  the 
notch  regioa  In  principle  the  stress  in  each  element  may  approach 
uniformity  by  decreasing  the  size  of  the  element.  This  will  increase  the 
number  of  the  elements  and  is  limited  only  by  the  computer  memory  and 
computation  time.  The  observation  of  symmetrical  model  is  obvious  from 
the  stress  values  of  full  specimen  model  and  half  specimen  model. 

B.      OBSERVATION  OF  LINEAR   RELATIONSHIP   BETWEEN  EXTERNAL   LOAD 
AND  LOCAL  TENSILE  STRESS 

The  tensile  stress  of  the  element  number  460  (Figure  4.3)  of  half 

specimen  model  which  subjected  to  the  displacement  B.C.  value  0.001  is 

0.002475  Msi  (Figure     4.9)  and  the  tensile  stress  of  the  same  element 

which  subjected  to  the  displacement  B.C.  value  0.004  is  0.0112  Msi  (Figure 

4.10).  The  ratio  of  tensile  stress  is  almost  the  same  as  the  ratio  of 

displacement,  the  different  value  bases  on  the  non-uniform  of  the  stress 

around  the  notch  region  and  can  be  accepted.  When  the  element  number  230 

(Figure     4.7)  of  the  same  model,  which  is  far  from  the  notched  region  the 

ratio  of  the  stresses  (Figure     4.11,  displacement  B.C.  value     0.001  and 

Figure    4.12,    displacement    B.C.    value      0.004)    and    the    ratio    of    the 

displacement  are  almost  the  same  and  the  difference  is  small  compare  to 

the  different  value  for  the  element  in  the  notch  region.  Therefore,  the 

linear  relationship  between  external  load  and  local  stress  can  be  accepted 
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whether  the  element  is  in  notch  region  or  not.  Since  the  element  is  small 
enough  to  be  satisfy  the  uniform  stress  within  element,  then  the  uniformity 
and  the  linearity  can  be  attained. 

C.      STATISTICAL  SIZE  EFFECT  RESULTS  DISCUSSION 

Equation  2.13  and  2.14,  calculated  the  reliability  of  the  specimen 
model  based  on  statistical  size  effect  in  two  dimensional,  non- 
homogeneous,  uniaxial  stress.  For  application  in  these  equation,  the 
displacement  boundary  condition  need  to  be  converted  to  the  average  stress 
along  the  boundary.  The  definition  of  the  average  stress  in  two  dimension 
is: 

(P}=  {  S  [wpqaxP<l] }/ {  S  wpq  } 

where:  pq  (element  number)  =  186,  187,  188,  189,  190,  200,  205,  210, 
215,  220,  225,  and  230  for  half  specimen  model  (Figure  4.7),  and 

pq  =  163,  164,  165,  166,  167,  168,  171,  174,  177,  180,  183,  186, 
189,  192,  465,  468,  471,  474,  477,  480,  483,  and  486  for  entire  specimen 
(Figure     4.5  and  Figure     4.13). 

The  values  of  stress  in  each  element  are  shown  in  Table  4.1  and  Table 
4.2    respectively. 

The  statistical  strength  of  a  structure  is  characterized  by  the  scale 

parameter   3^  in  equation     2.13  (  same  as  BETAA  in  computer  program, 

Appendix  A).  A  change  of  the  $E  value  is  therefore  related  to  the  strength 

change  of  the  structure.   The  values  of   0£  (equation    2.14)  for  plate  with 

notch,  where    the  values  of    <x,    £,  and    kc  are     10,    0.1  Msi,  and    0.001 
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0.001  respectively  are  shown  in  Table  4.  3-A.  Table  4.  3-A  also  shows 
that  the  variation  of  the  values  of  $£  which  are  calculated  from  different 
values  of  external  load  (PBAR)  from  different  values  of  displacement 
boundary  condition  (u').    The  values  of     3E  which  calculated  from  half 

specimen  model  seem  to  be  higher  than    0£  value  from  the  full  specimen 

model.  The  reason  is  that  the  number  of  elements  for  full  specimen  is  less 
than  the  number  of  elements  for  half  specimen  model,  therefore,  the 
accuracy  for  full  specimen  model  could  be  lower  than  that  for  half 
specimen  model.      For   the   purpose   of   reliability   calculation   for   the 

structure,  one  can  average  the    0£  values  from  one  set  of  finite  element 

data  output  from   different   values  of   load  under   the  same   boundary 

condition,  therefore.for  instance  the  average  value  of  £E  is  approximately 

0.030.  Table  4.  3-B  shows  the  values  of  0E,  where  the  values  of  <x,  and  kc 

are  changed  to  20,  and  0.0001  respectively.  The  0E  values  are  reduced  but 

the  difference  of  the  average     j$£  values  is  not  large,  because  of  the 

compensation  of    higher  <x  and    lower  kc  values.     This  means  that  the 

changes  of  <x  and  kc  values  do  not  affect  much  to  the  0£  value  but  <x  affects 
the  reliability  of  the  structure  (equation  2.13).  All  the  reliability 
calculation  results  for  different  values  of  £p  ,where  the  specimen  models 
are  operated  in  the  same  loading  range  are  shown  in  Appendix  B. 
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In  the  case  of  plate  without  notch  that  shown  in  Table    4.  3-C  and 

Table  4.  3-D,  the  values  of  £E  from  different  values  of  displacement 

boundary  condition  are  exactly  the  same.  The  reason  is  that  the  uniformity 

and  the  linearity  affect  the  accuracy  of  the  calculations.   The  value  of    $E 

from  the  plate  without  notch  (0.0439)  is  higher  than  from  the  plate  with 
notch  (0.028)  under  the  same  boundary  condition  (u'  =  0.01),  which  is  a 
reasonable  result. 

In  addition,  the  results  show  that  the  material  scale  parameter,  3  is 

converted  to  3E  by  the  geometry  of  the  structure  and  the  loading  boundary 

condition,  but  it  does  not  depend  on  magnitude  of  the  load.  The  results  of 
the  specimen  model  without  notch  confirmed  that  the  uniform  stress 
assumption  is  valid  and  known  result  can  be  recovered.  Between  results  of 

notched  specimen  model  we  showed  that  £E  remains  constant  for  different 

load  magnitude  because  the  uniformity  and  the  linearity  conditions  are 
attained.  As  expected  we  confirmed  the  notched  model  is  less  reliability 
than  the  plain  model  under  the  same  loading  condition.  Given  available 
experimental  data  of  a  composite  material,  we  can  evaluate  the  scale 
parameter  {J  and  shape  parameter  <x.  Using  these  results  we  demonstrated 
that  the  procedure  of  analyzing  the  effect  of  a  non-uniform  state  of  stress 
obtained  from  the  finite  element  method  can  be  post  processed  by 
statistical  reliability  formulation  accounted  for  size  effect,  resulting  on 
the  prediction  of  the  reliability  of  the  structure. 
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Figure  4.8  Stress  Distribution  Values  of  Half  Specimen 
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Figure  4. 13     Element  Numbers  of  Entire  Specimen  Model  at 
Bottom-Right  Boundary  Region 
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TABLE   4.1 
AVERAGE  5TRE55  FOR  HALF  SPECIMEN  MODEL  WITH  NOTCH 


ELEMENT 

ay  ,  Msi. 

a. 

NUMBER 

WIDTH,  IN. 

u'  =  o.oor 

U'  =  0.004" 

u'  =  o.or 

186 

0.057000 

0.009071 

0.036280 

0.090670 

187 

0.057000 

0.009358 

0.037420 

0.090350 

188 

0.057000 

0.009496 

0.037970 

0.094830 

189 

0.057000 

0.009554 

0.038190 

0.095370 

190 

0.057000 

0.009750 

0.038280 

0.095570 

195 

0.024875 

0.008758 

0.035030 

0.087550 

200 

0.024875 

0.008491 

0.033960 

0.084860 

205 

0.024875 

0.008179 

0.032710 

0.081700 

210 

0.024875 

0.007840 

0.031350 

0.078260 

215 

0.024875 

0.007500 

0.029980 

0.074810 

220 

0.024875 

0.007195 

0.028760 

0.071720 

225 

0.024875 

0.006968 

0.027840 

0.069400 

230 

0.024875 

0.006855 

0.027880 

0.068220 

PBAR* 

0.008717 

0.034850 

I 

0.087045 

Average  boundary  stress 
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TABLE  4.  2 
AVERAGE  STRESS  FOR  FULL  SPECIMEN  MODEL  WITH  NOTCH 


ELEMENT 

ay  ,  Msi. 

NUMBER 

WIDTH.  IN. 

U'  =  o.or 

163.  166 

164.  167 

165.  168 
171,  465 
174,  468 
177,  471 
180.  474 
183,  477 
186,  480 
189,  483 
192,  486 

0.095000 
0.095000 
0.095000 
0.024875 
0.024875 
0.024875 
0.024875 
0.024875 
0.024875 
0.024875 
0.024875 

0.091560 

0.094940 

0.095670 

0.087600 

0.084860 

0.081660 

0.078170 

0.074670 

0.071540 

0.069200 

0.068030 

PBAR* 

0.078025 

- 

Average  boundary  stress 
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TABLE    4.  3- A 
STATISTICAL  DATA  SUMMARY 

FOR  PLATE   WITH  NOTCH,  <x  =  10,   3  =  0.1,    kc=  0.001 


MODEL 

HALF 

FULL 

LIST       "^"-x^ 

U'=0.001" 

U'=0.004" 

u-o.or 

U'=0.01M 

PBAR* 

0.008717 

0.034850 

0.087045 

0.087025 

*E 

0.030582 

0.030749 

0.031110 

0.028880 

TABLE    4.  3-B 

STATISTICAL  DATA  SUMMARY 

FOR  PLATE  WITH  NOTCH,  ex  =  20,  3  =  0:1,    kc=  0.0001 


MODEL 

HALF 

FULL 

"^--^0iSPi.»csMer<iT 

U'=0.001" 

U'=0.004" 

u'=o.or 

u'=o.or 

PBAR* 

0.008717 

0.034850 

0.087045 

0.087025 

k 

0.029122 

_________ 

0.029298 

_ — _ __ __ 

0.029721 

0.028826 

Average  boundary  stress 
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TABLE    4.3-C 
STATISTICAL  DATA  SUMMARY 

FOR  PLATE  WITHOUT  NOTCH,  <x  =  10,  3=0.1,    kc=  0.001 


■ — — ^^DISPLACEMENT 
UlST.                  '     ~^ ~__ 

u'=o.oor 

U'=0.0l" 

PBAR* 

0.009512 

0.094910 

0.043948 

0.043948 

TABLE    4.  3-D 
STATISTICAL  DATA  SUMMARY 

FOR  PLATE  WITHOUT  NOTCH  ,  <x  =  20,  £  =  0.1,    kc=  0.0001 


*"~ — -^PfSPLACeMENT 

LIST                         "*— --^ 

u'=o.oor 

u'=o.or 

PBAR* 

0.009512 

0.09491 

k 

0.051800 

0.05180 

Average  boundary  stress 
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APPENDIX  A 
STATISTICAL  STRUCTURAL   RELIABILITY  PROGRAM 


C>I1IHII>>M|II    l>    MtOUl    MX>>|>IIM>4    IIKIIMIIMXIOIXI 

c: 

c 

c 

C  PROGRAM:  RELIABILITY  CALCULATION  OF  PLATE  MODEL  8ASED  ON 

C 

C  LINEARIZE  SIZE  EFFECT,  NONUNIFORM  UNIAXIAL  STRESS. 

C 

C  RELLA9ILITY  FORMULATION  IS  EXPRESSED  IN  THE 

C 

C  "    STANOAPO  tfELBULL    FORM    ". 

C 

C  DATA  INPUT:  OUTPUT  OF  FINITE  ELEMENT  CODE  NAME  "  NIKE20" 

C 

C»*****»*»**»**»*«»*»»«******«*»»***»*««»»*****»»*******«**»* 

c 
c 

C  VARIABLE:  KC  =  MATERIAL  CONSTANT  FOR  PLATE  THAT  IS 

C 

C  SUBJECTED  TO  UNIFORM  TENSILE  STRESS.  IN 

C 

C  THIS  CASE,  ASSUME  IN  EACH  ELEMENT  THE 

C 

C  LOCAL  TENSILE  STRESS  IS  UNIFORM. 

C 

C  BETA  =  SCALE  PARAMETER  FROM  EXPERIMENTAL  RESULT 

C 

C  FROM  SEVERAL  NuMRERS  Qe    SAMPLE. 

C 

C  ALPHA  =  SHAPE  PARAMETER  FROM  EXPERIMENTAL  RESULT 

C 

C  FROM    SEVERAL    NUMBERS    OF   -SAMPLE. 

C 

C  PBAR  =  AVERAGE  STRESS  ALONG  BOUDAPY  SIDE  OF  THE 

C 

C  PLATE  THAT  SUBJECTED  TO  ThF  EXTERNAL  FORCES. 

C 

C  IT  MUST  BE    CALCULATED  SEPARATELY  AFTER  OBTAINING 

C 

C  THE  STRESS  RESULTS  FROM  "NIKE2D"  PROGRAM. 

C 

C  N£  =  NUMBER  OF    ELEMENTS  OBTAINING  FROM  MESH  GENERATOR 

C 

C  PROGRAM  NAME  "MAZE"  WHICH  IS  THE  PROGRAM  THAT 

C 

C  CREATE  AN  INPUT  FILE  TO  BE  USED  IN  "NIKE2D" 

C 

C  PROGRAM.  ALSO  BOUNDARY  ELEMENT  NUMBERS  CAN  BE 

C 

C  CREATED  FROM  "MAZE"  PROGRAM. 

C 

c 

c 
c 

DOUBLE  PRECISION  F  ,  FA  ,FF  ,  P(,  lOOO>  ,  RA(  1  OOO)  ,  BETAA 
DIMENSION  XL(2363t  ,  XW(236a) .  S  I&U23fe8> 
REAL  KC 
BYTE  A<>30; 
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Tm*,  'Enter  file  containing  stresses  data  input • 
READC5,S)  A 
5   F0RHAT(30A1) 

0PEN<UNIT=*,  FILE=A»  STAnjS='OUD') 
OPEN(UNIT=6,  FIL£='SS.DAT',  STATIJS  = ' NEW ) 
\*RITE(6,1) 

1  FORMAT(3X. 'DATA  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT 
UNIAXIAL  TENSILE  STRESS') 

TYPE*, "ENTER   1   FOR  PLATE  WITH  NOTCH 

*  OTHERWISE  F0|'  PLATE  WITHOUT  NOTCH" 
READ*.NC 

IF(NC.EQ.l)  THEN 

«RITE(6.2) 
ELSE 

V*RITE(6,22) 
END  IF 

2  F0RMAT(2O*. 'FOR  PLATE  WITH  NOTCH') 

22  F0RHAT(18», 'FOP  PLATE  WITHOUT  NOTCH'! 

TYPE*,  'ENTER  THE  VALUE  OF  QISPLACENENT  BOUNDARY  CONDITION* 

PEAD*,U 

WRlTE(b,3)  U 

TYPE*, 'ENTER   1   FOR  ENTIRE  SPECIMEN  INPUT 

*  OTHERWISE  FOR  HACF  SPECIMEN  INPUT* 
READ*,NS 

IF(NS.EQ.l)  THEN 

WRI  I  £ ' 6-23) 
ELSE 

v»RITE(b»2f) 
END  IF  ' 

23  FORMATdSX, 'ENTIRE  SPECIMEN  MODEL  INPUT M 

24  FORMATf  16X, 'HALF  SPECIMEN  MODEL  IMPlfT'} 

3  FORMAT C7X. 'WHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  IS  *,F3.U) 
TYPE*, 'ENTER  THE  KC  VALUE' 
PEAO*,KC 

TYPE*»'ENTER  THE  BETA  VALUE* 
READ*, BETA 

TYPE*, 'ENTER    THE    ALPHA     VALUE' 
READ*, ALPHA 
ALPHA2=ALPHA+2.0 
FF=0.0- 

TYPE*., 'fNTER   THE    P3AR  VALUE' 
PEAD*,PBAR 

TYPE*,  'ENTER  THE  NUMBER  Op  ELEMENTS' 
PEAD*,NE 

C 

C      READ  ELEMENT  LENGTH,  WIDTH,  ANO' TENSILE  STRESS  FROM  INPUT  FILE 

C 

C      CALCULATE  BETAA  WHICH  IS  BETA  THAT  CHANGED  WITH  GEOMETRY  OF  THE 

C 

C      SPECIMEN  MODEL  AND  CERTAIN  LOADING  BOUNDARY  CONDITION  BUT  IT  LS 

C 

C       INDEPENDENT  OF  THE  VALUl  OF  LO^O  UNDER  TH£  SAME  BOUNDARY  CONOITION. 

C 

DO  200  1-1, NE 

READt4,20)     XLU),XWLI>.SIGYCI} 
IFCSrGY(I).GT.O.O)    THEN 

f  =  (xl(  rww'i)/Kc**2.o) 

»  *(S10T(I)**ALPHA2) 

FF=FP*-F 
ELSE 
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FF  =  FF 

END  IF 

2  00 

CONTINUE 

IFfNS.EQ.l)  THEN 
FF=PF 

EL36 

FF=FF*2.0 
END  IF 

FA=FF*«C-1 ./ALPHA2)»PBAR 
BETAA  =  8£TA**(  ALPHA/ ALPHA2)-*FA 
WRITE(b.lO)  KC 

10 

FORMATfax, 'KC  =',F\Z.9) 
WRITEC6.U)  &ETA 

1 1 

FORMA r ( &x. '8ETA  ='.Fi3.9) 
WRITE(6,i2)  ALPHA 

\? 

F0R«ar(5X, 'ALPHA  s',F!3.8) 
WRITE(6,13)  PBAR 

1 3 

FOW'Mf  f6x, 'PBAR  =',F13.8) 
WRI  rE«6, 1A)  NE 

i  a 

FORMAT C8K. 'NE  =  ' ,14) 
WRtTE(6,15)  ALPHA2 

IS 

FOR^ATCax, 'ALPHAS  :',F13.8] 

WRITE (6, 25)  BETAA 

TYPE*, 'ENTER  THE  FIffST  EXTERNAL  LOAD* 

R£aD»,P(1) 

c 
c 
c 

PIMC=Pf 1 )A30. 

PELTAftlLrTr  CALCULATION  BASED  ON  LINEARIZE.  SIZE  EFFECT 

DO  100  J=t,100 

J?A(J)=D£XP(.-<P<  J1/8ETAA^«*ALPHA2) 
WPIT£(6,30)  J,P(J),RACJ) 
P( Jtl )=P(J)*PJNC 
100   CONTINUE 
20   FORMAT(8X,2F11.0fE«2.0) 
25      FORMAT' 5X, ' BETA A  =' ,F\Z. 8) 

30   F0RMAP9X.  "T  =',I1,3X.'?    =•  ,F\1.Q,3X.  'RA    =',F13.3) 
sro  = 
END 
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APPENDIX  B 
DATA  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT  UNIAXIAL  TENSILE  STRESS 

DATA    OUTPUT    FOR    LINEARTZE    SIZE    EFFECT    UNIAXIAL    TENSll£    STRESS 
FOR    PLATE   WITH   NOTCH 
WHERE    Th£   DISPLACEMENT  BOUNDARY    CONDITION    IS         O.O01O 

half  specimen  model  input 


KC 

= 

0. 

OOOIOOOQ 

1 

8ETA 

s 

0. 

10000000 

1 

Alpha 

s 

20. 

OOOOOOOQ 

1 

P8AR 

s 

0. 

00871700 

N£ 

- 

48? 

ALPHA2 

: 

22. 

OOOOOOOQ 

1 

6ETAA 

~ 

0. 

0291221 2 

J 

s 

1 

P 

3 

'1.0087170  0 

RA 

3 

1 .00000000 

J 

= 

2 

P 

3 

0.10O10  757 

PA 

3 

1.000  0  0000 

J 

= 

3 

P 

3 

0.0  09  293  13 

RA 

3 

1 .00000000 

J 

s 

'1 

P 

3 

0.0095«"7<j 

Rrt 

3 

I  .  0  0  0  0  0  0  0  0 

J 

s 

5 

P 

3 

0.00987927 

RA 

3 

1.00  0  00000 

J 

= 

0 

P 

3 

0.0101b9«3 

RA 

3 

1.00000000 

J 

s 

7 

P 

3 

0.010460*40 

RA 

3 

1.00000000 

J 

3 

ft 

P 

3 

0.0  10  7509  7 

PA 

3 

1 .00000000 

J 

- 

9 

P 

3 

) .  0  1  1)  a  1 5  3 

RA 

3 

1.000000  0  0 

J 

z 

10 

p 

3 

0.0 1 1 3321 o 

RA 

3 

1 .00000000 

J 

3 

1  1 

p 

3 

0.01 162267 

RA 

3 

1 .00000000 

J 

z 

12 

P 

3 

0.01 1 91 32  3 

RA 

3 

1.00000000 

J 

S 

13 

P 

3 

0. 01220580 

RA 

3 

1.00000000 

J 

s 

ia 

P 

3 

0.0  12494  3  7 

RA 

3 

0.99999999 

J 

z 

15 

P 

= 

0. 01273493 

RA 

3 

0.99999999 

J 

= 

lb 

p 

3 

0.0130  7  550 

RA 

3 

0.99999O93 

J 

3 

1  7 

P 

3 

0.01 336607 

RA 

3 

0.99999996 

J 

s 

1H 

P 

3 

0.0  1 5  65  66  3 

RA 

3 

0.99999994 

i 

: 

19 

9 

3 

0.0 1394720 

RA 

3 

0.99999991 

.1 

z 

20 

P 

3 

0  .  0  1  4  2  3  7  7  7 

RA 

3 

0.Q9999OH5 

J 

s 

21 

P 

3 

0.01452333 

RA 

3 

0.99999977 

J 

s 

22 

P 

3 

0.0  14ft 1390 

RA 

3 

0.999999*5 

J 

3 

23 

P 

3 

0.O151O9U7 

RA 

3 

0.99999946 

J 

s 

2'x 

P 

3 

0.0154004  3 

RA 

3 

0.9999991  ft 

J 

z 

25 

P 

3 

0.0 15690*0 

pa 

3 

0.99999877 

J 

3 

26 

P 

= 

0.'H59M  1  7 

RA 

3 

0.99999ft  I  5 

J 

= 

?7 

p 

3 

0.01627173 

RA 

3 

0.99999725 

J 

3 

28 

P 

3 

0.01h5b230 

RA 

3 

0.9999959S 

J 

3 

29 

P 

3 

0.0!6353»»7 

RA 

3 

0.999994  0b 

J 

3 

30 

p 

3 

0.01  7 1  4  \  U J 

RA 

3 

0.999991 35 

J 

3 

31 

P 

3 

0.01  7  J  J 'JO" 

9  A 

3 

0.9999«747 

J 

3 

J2 

P 

3 

0.01  772  15  7 

RA 

3 

0.9999819ft 

J 

: 

33 

P 

3 

0  .  0  1  3  0  l  5  1  3 

RA 

3 

0.99997423 

J 

z 

3a 

p 

3 

0.0  1830570 

RA 

3 

0. 9999*336 

J 

3 

45 

P 

3 

0.01859627 

RA 

3 

0.99994819 

J 

3 

36 

P 

3 

0.018586*3 

RA 

3 

0.99992712 

J 

3 

37 

P 

3 

0  .  0  1  9  1  7  7  '4  0 

RA 

3 

0. 999ft9« 0  4 

J 

3 

ia 

P 

3 

0  .  0  1  9  '•  6  7  97 

RA 

3 

0.999«5ft0<3 

J 

3 

39 

P 

3 

0.01975^53 

RA 

3 

0.999R0337 

J 

3 

*0 

P 

3 

i).  02004  91  0 

RA 

3 

0.99972ft90 

J 

: 

Ul 

P 

= 

n..)20S3:>b7 

RA 

3 

0.99962797 

J 

3 

12 

u 

3 

0.02063  >  2  3 

RA 

3 

0.999  'i  9175 

J 

z 

«i 

P 

3 

0. 02092030 

RA 

3 

0.99R50870 

J 

3 

ua 

P 

3 

0.02121 1  57 

"A 

3 

0.99906375 

J 

3 

<4S 

P 

3 

0.0215  01<>3 

RA 

3 

0.99*73722 

J 

3 

46 

u 

3 

0.0P1 79250 

PA 

3 

0.9933037ft 

J 

3 

<*.7 

p 

3 

0  .0220830  7 

RA 

3 

0. 997730b! 

J 

3 

US 

p 

3 

^.(i2?573o5 

RA 

3 

0.99*97553 

J 

= 

U9 

13 

3 

0. 02266420 

RA 

3 

0.9959ft  n 76 
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DATA  OUTPUT  FOR  LINEARIZE  SITE  EFFECT  UNIAXIAL  TENSILE  STRESS 
FOR  PLATE  V«TH  NOTCH 
WHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  IS    0.0100 
HALF  SPECIMEN  MOOEL  INPUT 


KC 

= 

0. 

0  0010000 

BETA 

: 

0. 

10000000 

1 

ALPHA 

= 

20. 

oooooooo 

PBAR 

: 

0. 

O87O45O0 

1 

»£ 

= 

f82 

ALPHA2 

= 

22. 

00000 ooc 

1 

8ETAA 

3 

0. 

02972093 

J 

= 

1 

P 

3 

0.0  08  71  7  0  0 

PA 

3 

1  .00000000 

J 

S 

2 

P 

3 

O.i)  .1900757 

PA 

3 

1 .00000000 

J 

= 

3 

P 

3 

0.00<J?9<M  3 

RA 

3 

1  .  0  0  0  0  0  0  0  0 

J 

= 

4 

P 

3 

0.00958*70 

PA 

3 

1  .  0  0  0  0  0  0  0  0 

J 

z 

5 

P 

3 

9. 0  09 3 7 92 7 

PA 

3 

1  .00000000 

J 

= 

0 

P 

S 

0.01010383 

RA 

3 

1  .00000000 

J 

= 

7 

P 

3 

o.o t  ouooao 

RA 

3 

1.00000  0  00 

J 

S 

9 

P 

= 

0.01075)07 

RA 

3 

1  .00000000 

J 

3 

9 

P 

3 

O.oi 1 0U153 

RA 

3 

1  .00000000 

.1 

3 

10 

P 

3 

0  .  0  l  I  3  3  2  1  0 

RA 

3 

I  .00000000 

J 

3 

1  1 

p 

3 

0.011 622o7 

PA 

3 

1  .00000000 

J 

= 

12 

P 

3 

0.0  1 131 32  3 

RA 

3 

t .00000000 

J 

3 

13 

P 

3 

0.01220330 

RA 

3 

t .oooooooo 

J 

r 

1  <i 

P 

3 

0.012'!9'437 

RA 

3 

0.99999999 

J 

3 

15 

s 

3 

0.01278493 

RA 

3 

0.99999999 

J 

= 

lb 

P 

3 

0.01 307550 

RA 

3 

0.99999999 

J 

3 

17 

P 

3 

0.01 336607 

RA 

3 

0.93999998 

J 

s 

13 

P 

3 

0.01 365663 

PA 

3 

0.99999996 

,[ 

3 

19 

0 

3 

0.01  3TJ720 

PA 

3 

0 .93999994 

J 

- 

10 

9 

3 

0.01423777 

94 

3 

0.9999999  1 

J 

s 

21 

a 

3 

O.o  I  ■■152*3  3 

RA 

3 

0.99999985 

J 

= 

22 

0 

3 

0.0  1  as  I  890 

RA 

3 

0.93999978 

J 

3 

25 

P 

3 

0.01  5109'I7 

9  A 

3 

0.999999Ob 

J 

3 

?u 

P 

3 

0  .  i)  1  5  '1  0  0  (j  3 

RA 

3 

0 . 9999904" 

.1 

= 

?S 

P 

3 

0.0  1599060 

RA 

3 

0.99999921 

J 

= 

?6 

3 

3 

0.01538117 

RA 

3 

0.93999882 

J 

r 

?7 

P 

3 

O.i!  lb?  7  1  7  3 

PA 

3 

O.09999MPS 

J 

= 

2* 

P 

3 

0.0165&230 

RA 

3 

0.99999741 

J 

: 

?9 

P 

3 

0.ul6n5267 

RA 

3 

0.99999620 

J 

; 

30 

0 

3 

O.oi 714  5^3 

PA 

3 

0.93999447 

J 

= 

31 

P 

3 

0  .  0  1  7  '*  3  'X  0  n 

tf  A 

3 

0.99999200 

J 

s 

3? 

P 

3 

0.01 772157 

RA 

3 

0.99998848 

J 

s 

33 

P 

3 

0.01801513 

PA 

3 

0.99998553 

J 

s 

3" 

P 

3 

0.01 830570 

RA 

3 

0.99997658 

J 

s 

35 

P 

3 

0.01*59627 

RA 

3 

0.99996689 

J 

s 

3b 

P 

3 

0.01**8683 

RA 

3 

0.999O5343 

J 

: 

37 

P 

3 

1)  .  0  1  9  1  7  7  '4  0 

PA 

3 

0  .99993484 

J 

= 

38 

P 

3 

0  .  0  1  9  i  o  7  9  7 

RA 

3 

0.93990Q2', 

J 

s 

39 

P 

= 

0.019753S3 

RA 

3 

0.999*7434 

J 

3 

«0 

P 

3 

0  .  0  2  0  0  4  9  1  i.i 

RA 

3 

0.99982675 

J 

s 

ill 

o 

3 

0.  02033367 

RA 

3 

0.9997622'i 

J 

= 

ui 

p 

= 

0.020  6  302  3 

RA 

3 

0. 9996751* 

J 

= 

43 

? 

= 

0.02092080 

RA 

3 

0.99955817 

J 

S 

ua 

p 

3 

0  .  0  2 1 ? 1  13  7 

RA 

3 

0 .93940 1 5* 

J 

= 

*S 

p 

3 

0.021 50193 

RA 

3 

0.9991 92*4 

J 

3 

<46 

p 

= 

O.o?l  79JSO 

RA 

3 

0.9989157  0 

J 

3 

a7 

p 

3 

0  .  o  P  ? '!  rt  5  0  7 

RA 

3 

0  .9985491 h 

J 

3 

*8 

p 

3 

0.02237363 

RA 

3 

0.99806620 

.1 

3 

<*9 

p 

3 

0.02266'12*' 

RA 

3 

0.93743222 
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DATA  OUTPUT  F0«  LINEARIZE  SIZE  EFFECT  UNI 
FOR  PL ATE  WITH  NOTCH 
WHERE  TH6  DISPLACEMENT  BOUNDARY  CONDI 
HALF  SPECIMEN  MODEL  INPUT 


AXIAL  TEMSILE  STRESS 
TJON  IS    0.00*0 


KC 

= 

0. 

00010000 

l 

BETA 

s 

a. 

10000000 

ALPHA 

s 

20. 

oooooooo 

P8AR 

s 

0. 

03485000 

I 

N£ 

= 

482 

ALPMA2 

s 

22. 

oooooooo 

I 

8ETAA 

i 

0. 

02929771 

J 

s 

1 

P 

3 

9  .  U08  7 1  7 0 0 

PA 

3 

1.00000000 

J 

s 

2 

P 

3 

0.00900757 

PA 

3 

1 . 0  0  0  0  n  0  0  0 

J 

s 

3 

P 

3 

0.  J0OP93I 3 

9A 

3 

1.00000000 

J 

5 

n 

P 

3 

0.0O95  ^70 

9A 

3 

1 .oooooooo 

J 

5 

5 

P 

3 

0.0  09M7927 

PA 

3 

1.00000000 

J 

= 

6 

P 

3 

0  .  o  l  0  l  b  9  8  3 

9A 

3 

1.00000000 

J 

s 

7 

P 

3 

0  .  3  I  0 '•« 6 940 

PA 

3 

1  .oooooooo 

J 

s 

3 

P 

3 

0.0  I  0  7509 7 

PA 

3 

1.00000000 

J 

3 

9 

P 

3 

0 . 0  I  11  a  1 5  3 

RA 

3 

1.00000000 

J 

= 

10 

P 

3 

0.01  1  532.10 

RA 

3 

1  .oooooooo 

J 

= 

1  1 

P 

3 

3.01 1622b/ 

PA 

3 

1  .oooooooo 

J 

= 

12 

P 

3 

9. 01  191 323 

9A 

3 

t  .oooooooo 

J 

1 

)  3 

3 

3 

0.012293*" 

PA 

3 

1.00000000 

J 

3 

14- 

P 

3 

0. 1)1209137 

PA 

3 

0.9999999O 

J 

z 

15 

P 

3 

0  .  0  t  2  7  4  '4  9  3 

PA 

3 

0.99999OP9 

J 

= 

1  fc> 

u 

3 

O.Ot 307550 

PA 

3 

0,99909998 

J 

s 

1  7 

p 

3 

0.01  33*b07 

PA 

3 

0. 99999007 

J 

s 

18 

p 

3 

3.  01  3b5bb3 

PA 

3 

G.990PQ<aq«; 

J 

s 

19 

p 

3 

0.01  S9<172'J 

PA 

3 

O.PPPQPQOa 

J 

3 

20 

fc> 

3 

0..)  1  42  3  7  77 

Ob 

3 

0.999000*7 

J 

s 

?1 

p 

3 

0.01 452933 

PA 

3 

0.999999HO 

J 

s 

jj 

p 

3 

0.01  4  SI  «9" 

PA 

3 

0.9PP9PQb9 

J 

= 

25 

p 

3 

9.0  151 0  94  7 

PA 

3 

0. 90939953 

J 

r 

2" 

p 

3 

0.0  15  »00«3 

Pfl 

3 

0  .  9  9  9  0  9 '» 2  8 

I 

: 

?5 

p 

3 

O.'M  5 1°0 o-i 

PA 

3 

0.99900*0,? 

J 

= 

2b 

p 

3 

3.0159*1  i  / 

ft  A 

3 

0.9">999*3b 

J 

= 

?7 

p 

3 

J.0lo?717? 

PA 

3 

,).oqqqq75P 

J 

= 

?" 

p 

3 

9.0  1*5*230 

PA 

3 

0.9PP99b45 

.1 

= 

K. 

p 

3 

').J!6-*5207 

PA 

3 

0.qqqqq«79 

J 

r 

30 

p 

3 

3.01  7l'l  54  3 

PA 

3 

0  .99999P.42 

J 

= 

31 

p 

3 

0 . 0 1 7  '4  3  a  u  0 

PA 

3 

0  .93908003 

J 

3 

32 

p 

3 

0.0  1 7  72457 

RA 

3 

0.9O99H421 

J 

= 

33 

p 

3 

o.O 1801 51 3 

PA 

3 

0.OOpq7742 

J 

= 

m 

p 

3 

0. 91*30570 

PA 

3 

0.99pqb790 

J 

s 

35 

p 

3 

0.01859*27 

PA 

3 

0  .  qqqqvibU 

J 

s 

3fa 

p 

3 

o.i)18*H&93 

PA 

3 

0.99993615 

J 

= 

37 

p 

3 

0 .0 191 77uo 

°A 

3 

0.999OU)b7 

J 

z 

3  8 

p 

3 

0  .  0  1  a  (| ,_,  7  q  7 

PA 

3 

0.999H7564 

J 

3 

n 

p 

3 

0. 01975955 

RA 

3 

0.999*2772 

J 

3 

'4  0 

p 

3 

0.0  20  O't  9  1  o 

PA 

3 

O.9097b24« 

J 

3 

!»! 

p 

3 

9.0  29339*7 

Oi 

3 

0.999b7'«05 

J 

= 

42 

p 

3 

0.  02033^25 

RA 

3 

0.  99055'!  70 

J 

3 

'4  3 

p 

3 

0.  J20920MO 

PA 

3 

0.99039431 

J 

= 

'4-4 

p 

3 

0.02121137 

PA 

3 

0.9391 70*7 

J 

: 

'4  5 

p 

3 

0.O21501Q3 

PA 

3 

0.993?9357 

J 

3 

4h 

p 

3 

0.021 79250 

P4 

3 

0.99*51 375 

J 

3 

07 

p 

3 

9.i»22')8 3° 7 

PA 

3 

0.9980  1  1  '17 

J 

3 

aa 

p 

3 

0. )22373o^ 

PA 

3 

0.9973 '4  076 

J 

I 

»o 

p 

3 

0.022b*4Pi< 

on 

3 

0  .90«j4rt  1  SO 
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OATA   OUTPUT   FOR   LINEARIZE   SUE   EFFECT  UNIAXIAL 
FOR   PLATE   VtlTMOHT   NOTCH 
WHERE    THE    DISPLACEMENT    30UNOAR*.  CONDITION 
ENTIRE    SPECIMEN   MODEL    INPUT 


TCNSILF    STRESS 
IS         O.OOIO 


Kc 

s 

0. 

0001 oooo 

BETA 

= 

0. 

10900000 

ALPHA 

s 

20. 

oaoooooo 

P6AR 

z 

0. 

009S1200 

NE 

200 

ALPHAS 

z 

22. 

Oooooooo 

8ETAA 

3 

0. 

05180000 

J 

3 

1 

p 

3 

0.  00*71700 

pa 

3 

1  .00000000 

J 

S 

2 

p 

= 

0. .10900757 

RA 

3 

1  .00000000 

J 

3 

3 

p 

3 

O.OOopohi  3 

9ft 

= 

1  .  0  0  0  0  0  0  0  0 

J 

3 

U 

p 

3 

0.9095MH7O 

RA 

= 

1.00000000 

J 

3 

5 

p 

r 

0.')09M7927 

9A 

3 

1.00  0  00000 

J 

3 

fa 

p 

3 

0.0  toto9«3 

RA 

= 

1.00000000 

J 

3 

7 

p 

3 

o.oiouboao 

RA 

3 

1.0000  0  000 

J 

3 

« 

ra 

3 

0.010  75097 

9A 

3 

1.00000000 

J 

3 

9 

p 

3 

0.0  1  10015  3 

9A 

3 

1  .OOOOOOOO 

J 

3 

10 

p 

3 

0.011 33210 

RA 

Z 

1.00000000 

J 

3 

!  1 

p 

= 

0.01 Ib22b7 

9A 

Z 

1.00000000 

J 

s 

12 

p 

= 

0.01 191 323 

RA 

3 

1.00000000 

J 

3 

13 

9 

3 

0.01220380 

RA 

3 

1.00000000 

J 

3 

10 

p 

Z 

0.01249O57 

RA 

3 

1 .00000000 

J 

3 

15 

p 

Z 

0.  ;>  1  2  78  '19  3 

RA 

3 

1.00000000 

J 

= 

16 

p 

Z 

o.o I  30  7  550 

9A 

3 

1 .00000000 

J 

3 

17 

p 

z 

O.ili3ho07 

RA 

3 

I. oooooooo 

J 

= 

ia 

p 

= 

0.  0  1  3b5bh3 

RA 

3 

I  .  0  0  0  0  0  0  0  0 

J 

3 

19 

p 

3 

0  .  0  1  3  9  a  7  2  o 

RA 

3 

1.00000000 

J 

3 

20 

p 

= 

D.0  1  '423777 

PA 

3 

1. OOOOOOOO 

J 

3 

21 

p 

3 

0.0  1052 "3 3 

RA 

= 

1  .  0  0  0  0  0  0  0  0 

.  J 

3 

■>?_ 

p 

3 

,)  .  )  1  Oa  I  Hoo 

PA 

3 

1. oooooooo 

l 

3 

IS 

D 

3 

J.  i)  151  09*17 

RA 

3 

1  .  o  o  0  o  o  u  0  o 

J 

3 

2o 

p 

= 

0.  015 '•*  0003 

WA 

3 

1  .  0  0  0  0  0  0  0  0 

J 

= 

25 

p 

3 

il.!)l5b<»0bti 

PA 

3 

1  .oooooooo 

J 

r 

2b 

p 

3 

r».iil53ai  1  1 

PA 

= 

t .oooooooo 

J 

= 

27 

o 

3 

0  .  f>  1  o  2  7  1  7  3 

0\ 

3 

1.090000 0 0 

J 

3 

28 

? 

z 

u  .  i)  1  b  5  o  2  3  0 

RA 

= 

1  .OOOOOOOO 

J 

3 

29 

p 

z 

0."to^297 

PA 

= 

r.oooooooo 

J 

3 

30 

p 

z 

i).0  1  7  1  o  Vo3 

pa 

= 

1  .oooooooo 

J 

r 

31 

p 

z 

0.')  1  7  4  30  00 

WA 

3 

1  .  0  0  0  0  0  0  0  o 

J 

= 

32 

p 

= 

O.Dl  772'457 

i/A 

3 

1  .  0  0  0  0  0  o  o  o 

J 

3 

33 

p 

= 

0.01  HOI  51 3 

WA 

3 

1.000000 0  0 

J 

= 

S* 

p 

3 

0.0  1830570 

RA 

3 

1  .  00O00O0U 

J 

3 

35 

p 

3 

0.01H59b?7 

RA 

3 

1  .  0  0  0  0  0  0  0  0 

J 

3 

3« 

p 

3 

).  UH^S^ 

PA 

Z 

1.00000000 

J 

3 

37 

p 

= 

0  .01 o i  7  7  00 

9A 

Z 

I .oooooooo 

J 

= 

38 

p 

3 

0.01^0^797 

PA 

Z 

1 .oooooooo 

J 

= 

3» 

p 

= 

0.0 19 75 H5 5 

RA 

Z 

l .oooooooo 

J 

3 

<M> 

p 

= 

0.020  049  I ') 

RA 

Z 

1.00000  0  00 

J 

= 

0.1 

p 

3 

0.020  3  3967 

PA 

Z 

1 .oooooooo 

J 

3 

u-2 

p 

3 

0  .  0  2  o  *>  3  72  3 

PA 

Z 

1 .oooooooo 

J 

= 

a3 

p 

= 

0.02  0 52080 

RA 

Z 

l. ooo 00000 

J 

z 

aa 

p 

= 

0.021211  57 

PA 

Z 

1  .  0  0  0  0  0  0  0  0 

J 

= 

05 

p 

= 

0.021 50  1  >5 

WA 

Z 

1 .oooooooo 

J 

= 

a6 

p 

3 

0.0?1  7  9  2  ■a  0 

PA 

Z 

0'.99999ogq 

J 

= 

a7 

p 

3 

0.OPP0MJ07 

PA 

3 

0  .  ^popooqo 

.1 

= 

aS 

p 

= 

O.u22375o3 

PA 

= 

()  m  qqqqqoqq 

J 

= 

0.9 

u 

= 

0.!)2?hfi  12" 

RA 

3 

O.qqqqqqqo 

67 


OATA    OUTPUT    FOR    LINEARIZE    SIZE    EFFECT    UNIAXIAL   TENSILE    STRESS 
FOR   PLATE  WITHOUT    NOTCH 
WHERE  THE    DISPLACEMENT  BOONOABY    CONDITION  IS        0.0100 
ENTIRE    SPECIMEN   MODEL    INPUT 


KC 

s 

0. 

0OO1 0000 

SETA 

= 

0. 

10000000 

ALPHA 

: 

20. 

oooooooo 

PBAR 

; 

0. 

094-91000 

NE 

200 

ALPMA2 

: 

22. 

oooooooo 

8ETAA 

= 

0. 

05180001 

J 

= 

I 

p 

= 

0.0 08 71  7'}  15 

PA 

z 

1  .oooooooo 

J 

3 

2 

p 

3 

0.HO9 no  75 7 

PA 

z 

1  .  0  0  0  0  0  ft  0  0 

J 

S 

3 

p 

z 

0.. 10939°.  1  3 

PA 

z 

1  .  0  0  ft  0  0  0  0  0 

J 

3 

4 

p 

z 

0.')ii951H70 

PA 

3 

1.00000000 

I 

; 

5 

p 

z 

0  .009S792  7 

RA 

3 

1  .oooooooo 

J 

: 

b 

p 

z 

0.0  10  169*3 

PA 

3 

1  .oooooooo 

J 

= 

7 

p 

z 

0.0  1  046040 

PA 

3 

l.  oooooooo 

J 

S 

H 

0 

z 

0.  0  1 0  7509  7 

PA 

3 

1. oooooooo 

J 

; 

9 

p 

z 

0.01104153 

PA 

3 

I  .  0  0  0  0  0  0  0  ft 

J 

- 

in 

p 

z 

0.01135210 

PA 

3 

1  .oooooooo 

J 

s 

i  i 

p 

z 

0.01 tb2267 

PA 

3 

1.00000000 

J 

s 

12 

p 

z 

0.01 191323 

PA 

3 

1 .oooooooo 

J 

s 

13 

p 

z 

0.012203*" 

PA 

5 

1.00 0  0000" 

J 

s 

1  'l 

p 

z 

0.01239/1.57 

PA 

3 

1  .oooooooo 

J 

s 

15 

p 

z 

0.0l27»'i93 

PA 

3 

1  .oooooooo 

J 

z 

lb 

p 

3 

0.01  30  755"" 

PA 

3 

1.000000 0  0 

J 

s 

1  7 

p 

3 

o.OI 3  3os0  7 

PA 

3 

i.ooooooon 

J 

z 

18 

p 

3 

0.  )\ 3b5663 

PA 

3 

1 .oooooooo 

.!' 

s 

19 

p 

z 

O.OI 594 7P0 

PA 

3 

1 . o ft ooo oo 0 

J 

z 

20 

p 

z 

0.01123777 

PA 

3 

1  .oooooooo 

J 

z 

21 

p 

= 

0.01 452*3  3 

J  A 

3 

1  .  f,  0  0  0  0  0  0  0 

J 

z 

?2 

p 

z 

0  .  0  1  'HI  1 9  0 

PA 

3 

1  .oooooooo 

J 

z 

25 

p 

z 

0.0151 09a 7 

PA 

3 

1  .oooooooo 

J 

z 

?1 

p 

z 

0.01540005 

PA 

3 

l  .moooooft 

J 

z 

?s 

p 

z 

0.01  56906ft 

PA 

3 

1 .oooooooo 

J 

z 

•<» 

p 

z 

0.0159*1  t  / 

PA 

3 

1  .oooooooo 

J. 

z 

27 

p 

z 

ft.  iM62  7  I  73 

PA 

3 

1 .000  0  o  ft  0  0 

J 

z 

28 

p 

z 

0  .>  1  h5623  i 

PA 

3 

1 .oooooooo 

] 

z 

3Q 

? 

z 

0.016*52*7 

PA 

3 

1.0  0000000 

J 

z 

30 

p 

r 

ft  .11.71  a  3U  3 

PA 

3 

I  .  0  0  0  00  0  0  0 

J 

z 

i\ 

p 

= 

.i  .1 1  7-43'J  )') 

PA 

3 

I  .  0  0  0  0  0  0  0  o 

J 

z 

12 

k> 

= 

O.OI  772  <57 

PA 

3 

I. oooooooo 

J 

z 

33 

p 

3 

0.0  1*01 51 3 

PA 

3 

i. oooooooo 

J 

z 

34 

p 

= 

0.01830570 

PA 

3 

1.00000000 

J 

z 

35 

p 

3 

0.01959627 

PA 

3 

1. oooooooo 

J 

z 

56 

p 

= 

o  .0  l  «-sfls«  3 

PA 

3 

1  .oooooooo 

J 

z 

37 

p 

3 

0  .  ft  1  <»  1  7  7  4  0 

PA 

3 

1  .oooooooo 

J 

z 

3^ 

p 

= 

ft  .  'J  1  *» «6  1 9  > 

PA 

3 

1 .onoooofto 

J 

z 

39 

p 

3 

1.01975*5  3 

PA 

3 

1  .  0  3 0  0  0  0  0  0 

J 

z 

l»0 

c 

= 

0  .020049  1  0 

PA 

3 

1. oooooooo 

J 

z 

4-1 

p 

3 

ft.. IPO  3396/ 

PA 

3 

1  .OOOOOOOO 

J 

= 

»2 

p 

3 

0.020630  23 

PA 

3 

l  .oooooooo 

J 

z 

«3 

p 

= 

0  .  02r,:j2  18-) 

PA 

3 

1  .  oooooooo 

J 

z 

U4 

p 

3 

0  .  0  P 1 2  '  1  5  7 

PA 

3 

1 .oooooooo 

J 

z 

uS 

p 

Z 

0.021501  >3. 

PA 

3 

1 .oooooooo 

J 

z 

U6 

p 

z 

O.OPl  7935'' 

PA 

3 

0 .qqqnqoqn 

J 

s 

a? 

p 

r 

). <.i220*3.17 

PA 

3 

O.qsqagqoo 

J 

= 

a8 

p 

= 

n.  022 37 36 3 

P  A 

3 

0  .qpqPTJqq 

) 

= 

u>9 

p 

= 

O.o22b6'i2" 

"A 

3 

i]  „  399(3QT9'? 
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DATA   OUTPUT 


FOR   LINEARIZE   SIZE   EFFECT  UNIAXIAL   TENSILE   STRESS 
FOR  PLATE   with    NOTCH 
WHERE    THE    DISPLACEMENT    BOUNDARY    CONDITION    IS         0.0100 
ENTIRE    SPECIMEN   MODEL    INPUT 
0.00010000 


KC 

BETA 

ALPHA 

P8AR 

NE 

ALPHA2 

BETAA 

J 

J 

J 

J 

J 

J 

J 

J 

i 

J 
J 

J 
J 
J 

J 
J 
J 
J 
J 
J 
J 
1 

J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
I 
J 
J 
J 
J 
J 
J 
J 


=    3 


0.  10000000 
20.00000000 

0.03702500 
92 
22.00000060 

0.00280262 

I 

2 

3 

'4 

5 
fa 
7 

e 

9 
10 
11 
12 
1  5 
l'l 
15 
lb 
1  7 

18 

I  Q 
20 
21 
22 
25 
?u 
25 
2* 
27 
2« 
29 
30 
31 
32 
33 
31 
35 
3b 
37 
3« 
3° 
Un 

111 

U2 
u3 

(IS 
46 
a7 
ad 
a.9 
50 


P  = 

P  = 

P  s 

P  = 

P  = 

P  = 

P  = 

P  = 

P  = 

P  = 

P  = 

P  = 

B  = 

P  S 

P  = 

P  = 

P  = 

P  = 

P  s 

3  ; 

P  = 

D  = 

P  Z 


0.0087170  0 
3.0  090  0757 

0.0092911  3 
0.00958*170 
0.00987927 
0.01016983 
0.0  1  0  4<}040 
0.0  10  75097 
0.01 104153 
0.01 1 33210 
0.01 1622b7 
0.01191 523 
0.01220380 
0.0  124'-*"  37 
0.01?78'I93 
0.0  1 30  7550 
0.ni33tsb07 
0.01 3b5b»i 
O..I»1  394  72'1 
0  .  'J  1  a  2  3  7  7  7 
0.01  452853 
o.oi  'in  I  19.1 
0.0151 09«7 

0.015  4  0003 

0.0  I  569  "I  Mi 
0.015381  I  7 
0.0  I h271 73 
0.0  1  f>56230 
O.li  lb85.2H7 
0  .  0  1  7  1  4  5  0  3 
0.0  1 74  3000 
0.01 772457 
0.0180151  3 
0.0  I 8 505  70 
0.01 859627 
0.0  18°8b83 
O.O 191 7740 
0.01946797 
O.OI 975553 
0  .  0  2  0  0  4  9  1  u 
0.020 3396  7 
0.0 2 Ob  30 2 3 
0.02092T»0 
0.02121137 
0.02150  193 
0.021 79250 
0.  )2?'-*3"7 
0.0?2373b3 

0.  )22bb  *?• 

1 .  ii  2  2  ">  5  J  7  1 


RA 
RA 
RA 
PA 
RA 
PA 
PA 
RA 
RA 
PA 
PA 
PA 
PA 
RA 

HA 
PA 
RA 
C  1 
HA 
RA 
RA 
34 
RA 
<*A 
oa 

PA 
RA 
PA 
PA 
RA 
RA 
Pi 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
PA 
R  A 
RA 
RA 
R  A 
RA 
RA 

ka 

R  A 


I. 00000000 

1.00000000 
1.00000000 

1.00  0  00000 
I. 000 0  0000 

1 .00000000 

1 .00000000 
1 .00000000 

1 .oooooooo 

1.00000000 
1 .00000000 

1 . 0  0  0  0  0  0  0  0 

0. 99909990 

0.99909990 

0. 99999098 

0.99999997 
0.99999095 

0. 99909992 
0  .qqqqqqua 
0  .qqqqqqqi 
0.99999O71 

0 .93919954 

0  .gqqqqqiO 
0 . 9R909894 
0  ,999oq»t40 
0  .999M97bn 
0  ,99999b43 
g#qqgqq/l7U 

0  ,999oq22P 
0.9999887b 
0.99998373 
0.999O7fab0 
0.99996653 
0. 90995241 
0.99993271 
0. 9990053b 
0. 99906759 
0.99981567 
0.99974465 
0.99q(3u795 
0.93951688 
0.99934001 
0.99910233 
0.99878  4  27 
0.998  360 39 
0  .99779774 
1.99705383 
0.99607408 

'.'  .  O  3  ,'|  7  H  M  7  0 

0  . 9  9  J 1  0  8  8  1 
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DATA    OUTPUT    FOR    LTtfEARIZE    SIZE    EFFECT   UNI 
FOR  PLATE   WITH   NOTCH 
tfHERC  THt    DISPLACEMENT  BOUNDARY    CONDI 
HALF"    SPECIM6N    HOOEL    INPUT 


AXIAL    TENSILE    STRESS 
HON    IS         0.OO1O 


KC 

s 

0. 

ooiooooc 

9ETA 

s 

0. 

lOOOoooc 

ALPHA 

s 

to. 

oooaoooc 

I 

PBAR 

= 

0. 

00871 70C 

I 

N£ 

s 

482 

ACPHA2 

3 

12. 

oooooooo 

BETA* 

s 

0. 

03059189 

J 

3 

1 

P 

3 

0.  jr  B  7  1  70(1 

RA 

3 

0.99999971 

J 

s 

2 

0 

3 

0.000(10757 

PA 

3 

0.99999957 

J 

3 

3 

P 

= 

■j.  1)092991  3 

RA 

3 

0.99999°38 

J 

s 

a 

p 

3 

0.o09S«970 

RA 

3 

0.9399991 0 

J 

= 

5 

P 

3 

0.00997927 

PA 

3 

0.99999871 

J 

= 

6 

0 

3 

■0.01»lb993 

PA 

3 

0.9999991  7 

J 

s 

7 

p 

3 

0.0  1  r>'460ao 

(7  A 

3 

0.999997U« 

J 

: 

* 

p 

3 

0. .11075097 

Oft 

3 

0.999996Uq 

J 

s 

9 

p 

= 

1 .  0 1 1  0  a  ]  S  3 

RA 

3 

0.O9999S09 

J 

= 

1  ') 

s 

= 

0.011 33210 

RA 

3 

0.99999330 

J 

S 

1  1 

p 

= 

0.0  1  1  62  26  7 

RA 

3 

0.93999092 

.J 

= 

12 

p 

3 

0.0  1 131 32  3 

RA 

3 

0.99998779 

.J 

s 

1  5 

o 

3 

0.0  1220  390 

RA 

3 

0.93999369 

J 

- 

1« 

p 

3 

0  .  n  i  ?  -i  Q  'i  3  7 

RA 

3 

0.99997837 

J 

3 

15 

0 

z 

0.O127HUO3 

RA 

3 

0. 99997150 

J 

s 

1* 

p 

3 

0.01 307550 

PS 

3 

0. 99996268 

J 

z 

17 

p 

3 

O.oi33ob07 

OA 

3 

0.99995ia2 

J 

r 

1" 

p 

3 

0.01  3656o3 

RA 

3 

0.99993712 

I 

s 

19 

p 

3 

)  .  1  1  3  9 '4  7  20 

RA 

3 

0.9999190*4 

J 

s 

20 

p 

3 

0.0  1 U 23  7 7 7 

RA 

3 

0.999H9632 

.1 

s 

21 

p 

3 

0.01/152933 

PA 

3 

0.939867H7 

J 

s 

?2 

p 

3 

O.Ot  4"M  99" 

RA 

3 

0.999«32aa 

J 

s 

23 

(-• 

= 

0  .  0  1  5  I  0  3 0  7 

04 

3 

0.99978847 

J 

: 

?'J 

p 

3 

0.015 '4  0003 

RA 

3 

0.93973U15 

J 

3 

^5 

p 

= 

0.U15S906O 

RA 

3 

0.93o&<J732 

J 

= 

?6 

p 

= 

O.'MS^p  1  1  7 

w  a 

3 

0.93958539 

J 

= 

?7 

p 

3 

0.016271 73 

RA 

3 

0. 999085 3 a 

J 

s 

28 

D 

3 

O.'JI  6  5h2  3n 

pa 

3 

0.939}b359 

J 

3 

2g 

p 

3 

0.01o*52*7 

RA 

3 

0.99921595 

J 

= 

30 

p 

3 

o.oi7ia3a3 

RA 

3 

0.99903752 

J 

3 

31 

p 

3 

0  .  '"i  1  7  'i  3  '4  0  0 

34 

3 

0.93882257 

J 

~ 

32 

p 

3 

O.0  1  7  72 '45  7 

RA 

3 

0.9995o443 

J 

s 

33 

p 

3 

0.01901513 

RA 

3 

0.99825539 

J 

3 

in 

p 

3 

^.01830570 

PA 

3 

0  .99788649 

J 

s 

35 

p 

3 

0.0195992  7 

RA 

3 

0.9374 '4  740 

J 

s 

36 

□ 

3 

0  .  0  1 «  3  9  9  9  3 

RA 

3 

0.93o9262U 

J 

s 

37 

p 

3 

0.0W1  7  7'l" 

RA 

3 

0.996309 J6 

J 

3 

38 

p 

3 

0.019  1  fj  7  0  7 

RA 

3 

0.995581  10 

J 

= 

39 

p 

3 

0.01975*53 

RA 

3 

0.93472362 

J 

r 

ao 

a 

3 

0.02noa^io 

RA 

3 

0.99371655 

J 

= 

'4  1 

P 

3 

O.OPo  <;39fc.7 

!?A 

3 

0.99^53675 

J 

= 

u<> 

o 

3 

■1.0  20  9  302  3 

RA 

3 

0.931  1<S800 

J 

3 

u3 

p 

3 

u.  020 92 080 

34 

3 

0.98955067 

J 

3 

<ia 

-P 

3 

0.02124137 

*A 

3 

0.-99768135 

J 

= 

Q.5 

P 

3 

0.0?1501«3 

94 

3 

0.9*551251 

J 

: 

ab 

P 

3 

0.021 79  25 n 

RA 

3 

0.9*3  0  0210 

J 

= 

a7 

V 

3 

0.02209 30 7 

PA 

3 

0  .  9  8  0  1  0  3  1  8 

J 

s 

a* 

p 

3 

0.0?2373o3 

Ra 

3 

0.9767&3U9 

J 

: 

<*9 

r 

3 

O.0?26n  I2n 

w  a 

3 

0.9729251^ 
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DATA    OUTPUT   POR   LINEARIZE  SIZE    EFFECT   UNI 
FOR  PLATE    WITH  NOTCH 
WHERE  THE   DISPLACEMENT   BOUNOARY    CONOI 
HALF   SPECIMEN    MODEL    INPUT 
0.O01OOOOO 
.100OOOOO 
.00000000 
.03H850OO 


AXIAL   TENSILE   STRESS 
TIOH    IS  0.00*40 


KC 
8  ETA 

ALPHA 

P3AP 

NE 

ALPHA2 

8  ETA  A 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
.1 
J 
J 
J 
J 
J 
.J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 


0 

to. 

0. 
462 
12.00000000 
0.0307U947 

I 

2 
5 

4 

5 

6 

7 

8 

9 
10 
1  1 
12 
li 
1  i 
IS 
lb 
17 

18 

i? 

20 
21 

22 
23 
2< 
25 

?" 

37 

28 

29 
30 
31 
32 
ii 
3a 
IS 
Jo 
37 
38 
39 
U.0 
HI 
u2 
U3 
au 
IS 
ab 
1*7 
%6 

U9 


P  s 

P  = 

P  = 

P  s 

P  = 

P  s 

p  = 

O  ; 

P  S 

P  = 

P  = 

P  S 

P  s 

P  s 

p  = 

p  = 

P  s 

p  = 

p  = 

p  = 

p  = 

p  = 

p  = 

p  = 

p  = 

p  = 

p  = 

p  = 

p  = 

J  = 

p  = 

p  = 

p  = 

o  s 

p  = 

?  = 

p  = 

p  = 

p  = 

p  = 

p  = 

a  - 

P  = 

P  = 

P  = 

P  : 

o  - 

P  = 

P  : 


0  .  0  0  3  7  1  7  0  o 

0.0(>9y0757 

0. 0092981  3 

0.0  0'-»58  8  7  0 

0.009*7927 

0. 01016963 

0.0  1  O'4o0'40 

0.01075097 

0.0 1 1 oa 1 53 

0.01 1 33210 

0.01 1622b7 

0.01191 323 

0.  012203)30 

0.012'J9U37 

0.01273'i93 

0  .  0  1  3  0  7  5  5  i) 

0.01  336607 

0.01 365663 

0.01 39U720 

O.0  1U237  7  7 

0.01 aS2?33 

0.0  1  ai  I  890 

,0  15  I  0  9a 7 

,0  15  J 0  0 (i 3 

,01569060 

,015981  1  7 

,  n  |  <j  P  7  1  /  3 

,0165o23  ■) 

,016*529/ 

,ol  7i  'i  3  i  3 

,  0  1  7  1  3  a  0  0 

O.Di  772057 

0.0 1 6s t 51 5 

0.01330570 

0.  01859627 

0.0  !  «"-i633 

0.01917740 

0.0l°46797 

0.01975^53 

o .  o  2  o  o  a  9 1  o 

).02u33967 
O.'iPOs  502  3 
0.120  32  )«0 
l) .  0  2 1 2 1  l  5  7 
0  .  0  2  1  5  o  t  g  3 
0.021  7<*250 
0.02208  51!  / 
0.0  22 W 36 3 
D.02?-)6'I20 


PA 
94 
PA 
PA 
PA 
9A 
PA 
PA 

^4 

PA 
9A 
PA 

04 
PA 
PA 
PA 
PA 
PA 
PA 
PI 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
RA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
PA 
P4 
WJ 
PA 
PA 
PA 
PA 


0, 

,99999973 

0, 

,999999b0 

0, 

,999999^12 

I), 

,<?39q9qi5 

J. 

,99999«79 

0, 

,99999829 

0, 

,99999760 

0, 

,99999666 

0, 

.99999SU0 

0, 

,00999372 

0. 

.99999150 

0, 

,99998856 

0, 

,99993473 

o, 

,99997975 

0, 

,99997331 

0, 

.9999&505 

0, 

,99995450 

0, 

,9999ai 1 1 

0, 

.99992418 

0, 

.99990290 

0, 

.999*7626 

0 

.9998'4307 

0 

.99980 1 90 

0, 

.99975) 03 

0 

.99968843 

0 

.99961  1  70 

0 

.9P95  irtoo 

0 

,99Q'I0397 

0 

.99926570 

3 

.0990985* 

1) 

.99889726 

0 

,998655a9 

0 

.99336604 

0 

.99802.051 

0 

.99760923 

0 

.997121  Oh 

0 

.99654321 

0 

.99536101 

0 

.99505770 

0 

.  9  3  a  1  1  a  2  o 

0 

.99  311  08  8  0 

0 

.991  71 690 

0 

.990  2)067 

0 

.9*345373 

0 

.98ba2581 

0 

.98ao7233 

1) 

.9*1 35025 

0 

.97822221 

0 

.97462169 
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DATA    OUTPUT   PQR  LINEAR TZ£    SIZE    EFFECT   UN 
FOR   PLATE   WITH   NOTTH 
WHERE    THE    DISPLACEMENT 
HALF   SPECIMEN 


BOUNDARY    COND 
MOOEL    INPUT 


IAX1AL  TENSILE    STRESS 
ITION    IS         0.01OO 


<c 

s 

0. 

0O1OOO0C 

1 

SETA 

3 

0. 

1O0OOOO0 

1 

ALPHA 

3 

10. 

OOOOQOOfl 

1 

PBAR 

: 

0. 

0670460C 

1 

NE 

s 

*as 

ALPHA2 

s 

12. 

OQOOOOOO 

BETAA 

= 

0. 

03111025 

J 

= 

1 

P 

3 

0.00871700 

PA 

3 

0.99999977 

J 

s 

2 

P 

3 

0.00900757 

94 

3 

0.99999965 

J 

S 

3 

P 

3 

0.0092981  3 

R  A 

3 

0.99999949 

J 

= 

a 

P 

= 

0.00958870 

PA 

3 

0.99999927 

J 

s 

5 

P 

3 

0.00987927 

PA 

3 

0.99999895 

J 

s 

h 

P 

3 

0.0101o9«3 

RA 

3 

0.9999985  1 

J 

r 

7 

P 

3 

0.0  104604O 

PA 

3 

0.99999791 

J 

3 

8 

P 

3 

0.0t075097 

94 

3 

0.999O9710 

J 

= 

9 

P 

3 

0.0  1  104153 

PA 

3 

0.99999601 

J 

s 

10 

P 

S 

0.01 1 33210 

RA 

3 

0.99999454 

J 

: 

t  1 

P 

3 

0.01  162267 

PA 

3 

0.99999261 

J 

a 

12 

P 

3 

0.01 191 323 

RA 

3 

0.99999006 

J 

= 

13 

P 

3 

0.0122038O 

PA 

3 

0.99998672 

J 

s 

1  a 

P 

3 

0.01249437 

SA 

3 

0.99998239 

J 

S 

15 

P 

3 

0.01278493 

RA 

3 

0.99997680 

J 

s 

lo 

P 

3 

0.01  30755<> 

PA 

3 

0.99996962 

J 

; 

1  7 

P 

3 

0.01 3 36*0 7 

PA 

3 

0.99996045 

J 

; 

IJJ 

P 

3 

0.01  365bb3 

RA 

3 

0.99994880 

J 

S 

19 

P 

3 

0.01  394720 

RA 

3 

0.93995408 

J 

s 

?» 

P 

3 

0. OKI23777 

PS 

3 

0.999Q1558 

J 

= 

21 

P 

3 

0.O1  452833 

RA 

3 

0.999*9242 

J 

; 

22 

P 

3 

0.01481 890 

RA 

3 

0.99986557 

J 

= 

23 

P 

3 

0.O15I  0->47 

RA 

3 

0.99982777 

J 

S 

?a 

P 

3 

0.015*0003 

RA 

3 

0.99978354 

J 

3 

25 

P 

3 

0.0  15s90o'> 

PA 

3 

0.99972912 

J 

s 

?h 

p 

3 

0.01598117 

RA 

3 

0 .9996624 1 

J 

; 

n 

P 

3 

).')lb27173 

RA 

3 

0.99958  0  94 

J 

= 

;a 

P 

3 

0.01 65623^ 

RA 

3 

0.99948180 

J 

3 

29 

P 

3 

0.01685287 

RA 

3 

0.99936158 

J 

s 

30 

P 

3 

0.01714?  13 

PA 

3 

0.99O21627 

J 

3 

31 

P 

3 

0.0  1  7  4  3  '40  0 

RA 

3 

0.99904122 

J 

3 

32 

P 

3 

0.0  1772'45  7 

RA 

3 

0.998931 00 

J 

3 

33 

P 

3 

0.01301513 

RA 

3 

0.99857930 

J 

3 

3a 

P 

3 

0  .  0  I  8  3  0  5  7  0 

PA 

3 

0.99827883 

J 

; 

35 

P 

3 

0.01859627  ' 

RA 

3 

0.997921 17 

J 

= 

36 

P 

3 

0.0  1 388S83 

RA 

3 

0.997'I96ol 

J 

3 

37 

P 

3 

0.01917740 

RA 

3 

0.99699403 

J 

3 

36 

P 

3 

O.01<*46797 

RA 

3 

0.99640063 

J 

; 

3^ 

P 

3 

0.01975853 

RA 

3 

0.99570183 

J 

z 

ao 

P 

3 

0  .  0  2  0  n  a  -?  i  n 

RA 

3 

O.P9U88098 

J 

3 

III 

a 

3 

0.020*3967 

RA 

3 

0.99391915 

J 

; 

42 

P 

3 

0. 02063023 

RA 

3 

0.99279  4  8b 

J 

: 

<I3 

P 

= 

0.02092080 

RA 

3 

0.99148380 

J 

3 

(J  J 

p- 

= 

0.021 21 V3  7 

P4 

3 

0  *  90.995  p.55- 

J 

= 

US 

P 

3 

0.02150193 

PA 

3 

0.98818824 

I 

= 

46 

P 

3 

0.021 79250 

RA 

3 

0.9861 3821 

J 

3 

47 

P 

3 

0.02208  30  7 

RA 

3 

0.98376971 

J 

= 

as 

P 

3 

0. 022^7363 

RA 

3 

0  .9<M  0  594» 

J 

; 

tt9 

p 

3 

0.02266420 

Ru 

3 

0.97769°45 
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DATA   OUTPUT    FOR 


WHERE 

KC  = 

BETA  - 

ALPHA  = 

PBAR  = 

(IE  s 

ALPHA2  s 

SETAA  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  s 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 

J  = 


THE  DI 

EMTI 

0.00 

0.10 
10.00 

0.08 
592 
12.00 

0.02 

t 

2 

3 

a 

5 

6 

7 

a 
10 

1  1 
12 
13 
14 
IS 
16 
17 
IB 
19 
20 
?1 
22 
23 

2  'I 
2S 

26 
57 

PM 

39 

30 
31 

32 
33 
34 
35 
36 
37 
38 
59 

ao 

4,1 

aa 

43 
nu- 
15 
46 

y.  7 

U.M 


LINEARIZE    SIZE    EFFECT   UNI 

FOR    PLATE   WITH    NOTCH 
SPLACEMENT    BOUNDARY    CONOT 
RE    SPECIMEN    MOOEL     INPUT 
100000 
000000 
000000 
7025O0 


A*IAl    TENSILE   STRESS 
tion   rs       0.0100 


OOOOOO 
987978 

P  : 

P 
P 

P 

P 
P 
P 
P 

o 

p 
p 
p 
p 
p 
p 
p 
p 

p 

p 

o 

p 

p 
p 

p 
o 

p 
p 
p 
p 
p 
p 
p 
p 
p 
p 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 

o. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 

0. 

0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 


00871 700 
00900757 
00939*1 3 
0O953870 
009*7927 
01016983 

0  1  0  4  &  J  4  0 
01075097 
01 104153 

01  133210 
01 1622o7 
01 191 323 
01220330 
01249437 
01273  49  3 
01 307550 
01  336607 
OI36S6o3 
0  13  9  4  7  20 
01423777 
01452*33 
01431*90 
0  151 094  7 
0  1  5  '» 0  0  0  t, 
0156906''* 
0159H1 1 7 
016271  73 

0  1  g5o2"^0 
01635237 

01  714343 
0  17  4  3  4  0  0 
01772457 
0180151 3 
01630570 
01359627 
0  1  8*85*3 
01917740 
U  t  9  '1 6  7  9  7 
01975*53 
0  20  0  4  9  10 
020*3967 
0  2)63021 
0  20  9  2  0  3  0 
02121 137 
021 50 193 
021  7925't 
0220330  7 
02237563 
"226o  J2o 


94 
PA. 
RA 
PA 
RA 
PA 
RA 
PA 
RA 
PA 
RA 
RA 
RA 
RA 
RA 
RA 
Oft 
RA 
RA 
RA 
RA 
Oft 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
PA 
RA 
RA 
RA 
PA 
RA 
RA 
RA 
RA 
04 
PA 
RA 
RA 
■5  A 


0.99999943 
0.99999915 

0.99999376 
0.99999321 
0.99999743 
0.99999636 
0.99999490 
0.99999292 
0.99999025 
0. 99998608 
0.99998195 
0.99997572 
0.99996753 
0.99995700 
0.99994334 
0.999O2581 
0.99990342 
0.99987498 
0.99983905 
0.99979J87 
0.99973733 
0.999be>o89 
0.999579/19 

0.99947153 

0.99935869 
n.9991 7567 
0.9939770  4 
0.9997J51 ^ 
0.99944 1 9 1 
0.99903737 
0.997660U7 
0.99714793 
0.99653447 
0.99530P45 
0.99495150 
0.99389825 
0. 99267590 
0.99123333 
0.9*953721 
0.99754645 
0.°*521678 
0.992  4  9774 
0.97933261 
0.97565798 
0.97140315 
0.96648977 
0.960831 39 
0.95453317 
0.94639176 
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DATA  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT  UNIAXIAL  TENSILE  STRESS 
FOR  PLATE  WITHOUT  NOTCH 

WHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  IS    0.001O 

ENTIRE  SPECIMEN  MODEL  INPUT 

ICC  =  0.OO1  00000 

BETA  =  0.10000000 

ALPHA  =  10.00000000 

pBAR  =  0.00951200 

ME  =  200 

ALPHA2  =  12.00000000 

BETAA  =  0.0*394803 

J=  1  P  =    0.00871 7)0  R4  =  1.00000000 

J  =  2  P  =    0.0o90j757  RA  =  0 .  99999999 

J  =  3  p>  =    0.)0924?13  "J  =  0.9999999O 

J  =  a  °  =    C. 00958870  °A  =  0.9^909490 

J  =  5  P  =    0.00987927  RA  s  0. 93999098 

J  s  6  P  =    0.1t01o983  "4  =  0. 99999098 

J  =  7  P  =    0.0  1  U'ibO'lO  96  =  0.99999^97 

J  =  8  P  =    0.01075)97  94  =  0.93999995 

J  =  9  P  =     0.D1  t  0  a  I S  3  9A  =  0. 9999949a 

J  =  1 0  P  =    0.01133210  RA  =  0.999999O1 

J  s  11  P  =    0.11162267  R4  =  0.93999988 

J  =  12  P  =    0.01191323  RA  =  0.99999o<3U 

J  =  13  P  =    0.01220381  RA  s  0.99999979 

J  s  1 ')  o    -          0.112'»9'i37  PA  =  0.99999972 

J  =  15  P  =    0. 01278493  RA  =  0.93994963 

J  =  lo  t>  =    0.0130755"  RA  =  0.939O9952 

J  =  17  ^  =    0. 01336607  R4  =  O.99999037 

J  =  18  P  =    0.013b5s63  RA  =  0.99999O|0 

1=13  P  =    0.01394720  R4  =  0.99909396 

J  =  20  P  =     0. 1)023777  h>A  =  0.93999*66 

J  =  21  P  =    3.01'lS2955  RA  =  0.99999830 

J  s  2  2  P  =    9.Q1<i3l<t9<)  94=  0.999  99  78a 

J  =  23  °    =    0. 01510947  94  =  (j.93999727 

J  s  >'4  o  =    O.OlSaCOOJ  04  =  0.99999657 

J  =  25  P  =     (.0153906')  RA  =  0.999O9571 

J  =  2*5  P  =    0.0159M17  PA  =  0.99990465 

I  =  27  P  s    0.01627173  RA  =  0.999O9336 

J  =  2a  P  *    0.01656230  pa  =  0.9999917O 

J  =29  ,J  =    0. 010^52^7  94  =  0.99998989 

J  =  30  0  =    0.')1  7  1  'I  343  ^A  =  0. 99998750 

J  =  31  p  =    0.0  1  7'13'iOO  94  =  0.99998481 

J  =  32  P  =    0.01772457  RA  =  0.°999B148 

J  =  33  P  =    0. 01801514  °A  =  0.99997749 

J  =  34  p  =    0.01330571  RA  =  0.93997273 

J  =  35  P  =    0.niP59«s27  R4  =  0.99996705 

J  =  36  o  =    0.01888683  94  =  0.93996031 

J  s   37  P  s    0.01917740  94  =  0.9999523a 

J  =  38  P  =    0.019-46797  94  =  0.999Oa201 

J  =  59  P  =    0.01975*53  RA  =  0.99993180 

J=  4  0  P=     0.0 200 4 910  94=  0.9  999  18  75 

J  =  ai  P  =    0.02033967  3A  =  0. 09990343 

J  =  42  p  =    0.02065023  94  =  o.999«8551 

J  =  43  P  =    0.0?09P08'1  94  =  0.909*6459 

J  =  UU  P  =     0.03|>!1-37  94  s  !).999«4022 

J  =  45  P  S     0.02150143  R4  =  0.999KM89 

J  =  46  P  =     0.ii?1792S0  34  =  0. 99977901 

J  =  U7  p  =    i».0?20«307  R«  =  0.O3974095 

J  =  48  P  =     1.02237365  =4  =  0.93969646 

J  =  49  P  s    0.0?26o'420  PA  =  0.°9964622 
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Data  output  for  LINEARIZE  s*ze  effect  UNI 

FOR   PLATE    WITHOUT   NOTCH 
tfHERE    THE    DISPLACEMENT    BOUNDARY    CONDI 
ENTIRE    SPECIMEN   MODEL    INPUT 

o.ooiooooo 

0.10000000 

10.00000000 

0.09491000 
200 
12.90000000 

0.04390302 


AXIAL  TENSILE    STRESS 
HON   IS         0.0100 


KC 

BETA 

ALPHA 

P8AR 

NE 

ALPHA2 

8ETAA 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

J 

I 

J 
J 
I 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
J 
I 
J 
J 
J 
J 
J 
J 
I 
J 


l 

2 
3 

q 

5 

6 

7 

8 

9 

1  0 

1  I 

12 

13 

1  4 

15 
lo 
1  7 
l« 
|9 
?0 
21 
22 
25 

?y 

?5 
26 

27 
?9 
29 

3" 
31 
32 
33 
3 'i 
35 
3h 
37 

3° 
ao 
ai 

u3 
au 
'J  5 
(id 
U7 
48 
49 


P  = 

D  - 

P  = 

P  = 

P  = 

P  = 

P  = 

P  = 

P  = 

P  E 

P  = 

P  = 

P  = 

P  S 

P  = 

P  = 

P  : 

p  = 

P  = 

p  = 

P  ; 

P  = 

P  = 


0.00871  701) 
0.00900757 

0..)i)9SH"*7  0 

o.  inq°7'>27 
0.0101h9«3 
0  .  0 1  0  '4  6  0  4  0 
"0.U1  075097 

o .  o  1 1  i  a  t  s  5 

0.11  1  33210 
0.01  162267 
0.01  191 323 
0.  0122038  U 
0.012  4943  7. 
0.0127 -3  a  93 
0.01 307550 
O.'Jl  336607 
0.01  3656e3 
0.0  1  3  '  '1  7  ?  0 
j.Oia-2  3777 
0.01 '452933 
0.0 1  -IS  1*90 
1.0  151094  7 
0.01540003 
0  .  ii  1  5  b  9  Toi 
0.0  15961  17 
0.016271 73 
0.0  1  o5o23o 
0.01o?52"7 
> .  a  i  7 1  n  3  n  5 

0.01  7  4  3 '» 1  ■'" 
0.01  772457 
0.01801513 
0.O1A3057O 
'>.0l85-9b27 
0.0  1**«6-»3 
0  .  0  1 9 1  7  74  i 
0.0  1« '4  6797 
,  0.1-9  75953 
,  0  ?  0  0  i  9  1  u 
,02u339o7 
,0?<>&302  3 
,  r):?O92'M0 
,i21 ?1  1  37 
,.121^.1 1  '->$ 
•).02l  79^5:i 
O.0?20*30  7 
"1.022  573n? 
).'i2?'j;viPO 


PA 
94 
9  4 
9ft 
RA 
RA 
PA 
PA 
RA 
Rft 
RA 
PA 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
RA 
PA 
PI 
PA 
RA 
RA 
RA 
PA 
Rft 
RA 
"(A 
RA 
RA 
RA 
PA 
R  A 
9A 
PA 
PA 
PA 
RA 
RA 
PA 
RA 
^ 
P  A 
PA 
^.\ 
RA 
p  i 


1 . 

oooooooo 

0, 

qqqqqqqq 

0. 

gqgogqqQ 

0. 

qpgqgeigq 

0. 

qqgqgqgH 

0. 

qqqqqqqp, 

0. 

,qgqqq997 

0. 

^qqo^'iqS 

0. 

tqqqqqqqU 

0. 

,qqq.iqqqi 

0, 

,qoqq9988 

0. 

qqqqqq^H 

0. 

qqqqqqjq 

0, 

,99999972 

0, 

,99999963 

0, 

tqqqqqq=,^ 

0, 

(gqgqq9i7 

0, 

tqqqqqq  1  9 

0, 

.999998^6 

0, 

,9P999Hbh 

0, 

.99999850 

0, 

.99999794 

0. 

.99999727 

0, 

,qqqqqfc57 

0, 

,gqq.)gS7| 

0, 

.9999qu05 

0 

.99999336 

0 

.999991 79 

0 

.99998919 

ll 

,qqq0H759 

0 

.  ?9998481 

0 

.99998148 

0 

.99997749 

0 

.99997273 

0 

.99996705 

0 

.99996031 

0 

.99995234 

0 

,<3qqqnpq\ 

f) 

,q3qq3l HO 

0 

.99991875 

n 

.99990343 

0 

.9991R551 

0 

.99986459 

0 

,99994022 

0 

.9991 1 1 89 

0 

.9  99  7  790 1 

1 

.99974095 

0 

.903*,qh9(, 

0 

.994,4622 

75 


APPENDIX  C 

COMPUTER  PROGRAMS  FOR  CONVERTING 
FINITE  ELEMENT  (NIKE2D)  PROGRAM  OUTPUT  TO  BE 

STATISTICAL  RELIABILITY  COMPUTER  PROGRAM  INPUT 

c 

c 

C      PROGRAM:  CONVERT  NIKE2D  OUTPUT  TO  BE  STAT1S1CAL  INPUT 

C 

C  FOR   1PT.    INTEGRATION    ENTIRE    SPECIMEN  /10O6L   ONL* 

C 

c***»»*»«**»*»**»»»»»»»*»«»****»»»*»*»**»*»»*»«»»»*»»***«**»«* 
c 

c 

DIMENSION    Xl_(59<u  ,XW(5940. 
3YTE    A(30-) 

TYPE*, 'ENTER    FILE    CONTAINING    NIKE2D    DATA    OUTPUT' 
READ    (5,500)   A 
500    FOHMAT/30AI) 

0PEN(UNIT=S.     FILE=A,    S TATUS=  'OUO  '  1 
0PEN(UNIT=6,    FILE='CC.DAT' ,     STATUS= ■ NEW ' ) 
HR1TE(6,S0) 
NE=S92 
NP=NE/10+1 
00    100    M=l  ,  18 

*L(M)=. 1658354 
XW(M)=.0  95 

100  CONTINUE 

DO  101  M=19,36 
XL(M)=.0995 
XW(M)  =  .0  9S 

101  CONTINUE 

00  10?  M  =  37,  132 
XL(M)=. 024875 
Xrf(M)=.095 
10?   CONTINUE 

00  1  03  M=133,150 
XLCM)=.0995 
XW(M>=.095 
103   CONTINUE 

00    104  M=151 ,168 
XL  (M)  =  . 165833 <+ 
XW(M)=.09S 
10a       CONTINUE 

DO    105    M=169,200 
XL  CM)  =  .1658334. 
XW(M)=.02487S 
105      CONTINUE 

DO    106   M=20i ,216 
XL(M)  =  .09<J5 
XW(M)  =  .02«f875 
10b       CONTINUE 

DO  107  M=217,232 
XL (Ml =.024875 
X*(M) =.024875 

107  CONTINUE 

DO    1 08    M  =  233,296 
XL(M>=.099S 
XW(M) =  .02<f875 

108  CONTINUE 

DO  l  09  f*-?97  ,  316 
XL  CM)  =  . !  653334- 
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<W<M)=.02Ma7S 
I'll   CONTINUE 

DO  110  (V»=317,328 
XL  CM)  =  .1658334 
XW(M)=.02487S 

110  CONTINUE 

00    111    N=329,360 
XtfMl=.099S 
X*(M)=. QZH&7S 

111  CONTINUE 

DO  J  12  N  =  J61,42U 
XL  (M)  =  .02<*875 
XW(M>=. 024875 

112  CONTINUE 

00    113   M=a2S,456 
XL(N)s.0995 
XW(Mj =.02*875 

113  CONTINUE 

DO   114  M=457,438 
xl (M>=. 1656334 
X*(M)  =  .024875 
11 'l      CONTlNue 

00    1  IS   M=<f&«»,576 
XL (M)=. 024875 
XW(Mi =.024875 
115      CONTINUE 

DO    1  16   H=577,592 
XLfM)=. 02*875 
X*(M>  =  .02487.5 
lib      CONTINUE 

00    10    i  =  i  ,NP 
DO   15    7=1  ,10 

M=10*(I«!)+j" 
IFlW.Nt.NE'    7H£N 

JFfJ-.Eq.10)    THPN 

Nt"AD(5»45)     SIGY,SIGZ,SI3X.SI«YZ 
ELSE 

REA0(5,35)    SIGY*5TGZ,SISX,SIGYZ 
END     IF 
ELSE 

READ<.5.30)     3IGY,SIGZ,SIGX,SIGYZ 
GO    TQ    55 
END    IF 

WR[TE,(6.40)    M»XL(M)#XH<«) ,  S IGY ,  SI  GZ,,  S IGX ,  SIGYZ 
15  CONTINUE 

1  <j    CONTINUE 

SS    WRITE(e>,«fO)    M,XLCM)  ,XW(Mi  ,SICY,SlGZ,SIGx,SIGYZ 
jo    FORMAT( i8X,4£U .0) 
55   FORNATf 18X,4E11 .0.////) 
4')    FORMAT (,3X, IS, 2F 1 1.7,4*12-4) 
<*5    FORMAT  f  18< ,4E1 1  .0,  /////// /////) 

50    FORMAT C IX.  'ELEMENT   NO. • , »X, • L ' , 1 0X, "W ' ,7X, ' SIGY • ,8X , ' SI GZ 
.,8X, 'STGX' ,7X. 'SIGYZ' ,/) 
STOP 
END 
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c *....»................».. ... ........... .......... .......... 

c 

c 

C      PROSRAH:  CONVERT  NIKE20  OUTPUT  TO  B£  STATISTICAL  INPUT 

C 

C  FOR  4PT.  INTEGRATION  HALF  SPECIMEN  MODEL  ONCY 

C 

c*»**.*»*»*»*»».**«.***»»»«»*»*. *»»»»«*»»»**»««*»*#*»»»*»«#* 

c 

c*«»»»*»*»*«*. •««»«*». »»**•»»>«****•**«»*»**•►****»»»«»*»#•* 

c 

DIMENSION    XL('9?tf)  ,yvt(\92S) 
BYTE    A(30) 

TyPE*,'£NT£R    rILE    CONTAINING    NIKE2D    OUTPUT' 
REAO(5,500)    A 
500    FORMATC50A1  ) 

OPEN(UNir=S,    PlLE=A»     STaruS='OLt>'  ) 
OPEN(UNTT=6,     FILE=' W.Oar • ,    S F A TUS= *  NEW' ) 
WPt TE(b,50) 
NE=*8? 
XE«=ME«4 
NP  =  Ne4/40*1 
00    1 OO    M=l,100 
XL(H)  =  .0<t975 
XW(M)=. 02850) 
100       CONTINUE 

DO   iOi   M=101,220 
XL  ff«)  =  .  024875 
XW<M>=.O2Q50O 
10  1        CONTINUE 

DO   102  "=21 1»540 
XU*n  =  . 0124375 
XW<f1i  =  .0295  000 

102  CONTINUE 

DO  1<J3  M=541  ,660 
XL (*)  =  . 02^875 
XWfM)=. 028500 

103  CONTINUE 

00    i  0^  Mi  =  661,760 
X I (p)  =  .04975 
xW(M)=. 02850 

104  CONTINUE 

DO    105    M=761 ,920 
XLfM)=. 04975 
XW(M) =.0124375 

105  CONTtMUE 

00     lOfc  M  =  921  ,1112 
xL(M)=.02afl75 
XW(f4)=.0l2'J3  75 
10b      CON'l^'UE 

00    107  M=l 1 13, 1368 
XL(M)=. 0124375 
XWfM)=. 0124375 

107  CON  F  1WE 

00  1»«H=  1369,  1560 
XLfM)=. 024875 
XW(M)=. 0124375 

108  CONTINUE 

DO   tO<}   M=I5«),>720 
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XL (M> =.04975 
XW(W)=.012437S 

109  CONTINUE 

00    110    M=1721 ,1896 
XL(M1=. 01243 75 
XW(M)=.0l2437S 

110  CONTINUE 

DO  111  M=1897,t928 
XL (M) =.0124375 
XW(M)=. 0124375 

111  CONTlMUt 

00  10  I=1,NP 

00  15  J=1 » I o 

00  20  <=  1 ,4- 

MsftO»(I-1  )>4.*(J-n+K 
IF(*.NE.N£4)  THEN 
IF(K.EQ.«)  THEN 

TF(J.EQ.IQ)  THEM 

READ(S,4.S1  SIGY,3IG7:,SIGX.,5IGYZ 
ELSE 

READ(5,3S)  SIG*,SIGZ,SI6X,SIGYZ 
EMO  IF 
ELSE 

READ(5r30)  SlGf,SlGZ.SIG*,SlGYZ 
ENO  IF 
ELSE 

READ'S, 30)  SIGY,SIGZ,STGX,STGY7 
GO  TO  55 
ENO  IF 

WRITE (6,40)  M,XLCM) . XWCM) . SIGT ,  STGZ, SIGX, S IGYZ 
20        CONTINUE 
IS     CONTINUE 
10  CONTINUE 

55  WRITEf6,40)  M,XLCM},XWCM) ,3!Gy,5IGZ,SIGX,SIGYZ 
.30  F0RMAT(iax,uE11  .0) 
35  FORWAT{ 18X,4E11.0./) 
40  FORMAT  (3x.,T5f2P11  •  7-4E12.4) 
us  FQR.MAT(t9X,4E-H.0,/////////J 

50    FORMATftX,  'ELEMENT  NO . '  , 4X,  • L ' ,  1  OX,  ■ W  ,7X,  ' S IG Y ' , dX ,  ' S  IGZ ' 
*,3x, 'SIGH' ,7X, 'SIGYZ* ,/) 
3T0P 
ENO 
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c 
c 

C  PROGRAM:    CONVERT  NTKE2D   OUTPUT  TO   BE   STATISTICAL    INPUT 

C 

C  POR   iPT.    INTEGRATION  HALF  SPECIMEN  MODEL    ONLt 

C 

c 

c 

DIMENSION    XH4821  »XW(402) 
B>T£    A(iO) 

TyPE*,'FIL£    CONTAINING   NIKE20   DATA    OUTPUT' 
READ(S,5O0)    A 
500    FORMATC30AD 

OPEN((jNtT=S,    FIUE=A»    STATuSs'OLO') 
0PEN(UNIT=6,    FlLEs'aa.OAT' »    STATUS* 'NEW  ) 
WRITE(b,50) 
NE  =  4-82 
NP=*6/10*1 
00    100    Msi ,25 

XLfMl=.099S 

*"y»(m)  =  .0570 

100  CONTINUE 

DO    10 1    M=2$,SS 
XUM)=. 04975 
XW(M)=.0S7 

101  CONTINUE 

00    102  M=56,13S 
XL.fM)  =  .  02*875 
XwfMi=.4S70 

102  CONTINUE 

00    103   M=136rt65 
XHM):.  04975 
XWM)  =  .0570 

103  CONTINUE 

00  100.  M=t66,  190 
XLfM)=.099S 
XW(M)=.057O 
10U   CONTINUE 

DO  105  M= 191 ,230 
XL(M)=.09*5 
XV»(V»)  =  .024875 

105  CONTINUE 

DO    10b    M=231 ,278 
XLCN)s. 04975 
XW(M)=. 024875 
10b       CONTINUE 

00    107    M=279,3o2 
XL (M) =.02*875 
XW(M)=. 024875 
107       CONTINUE 

00    10*    M=3<K3,390 

XLfM)  =  .0^-97S 
X\N<M)  =  .02«*875 

106  CONTINUE 

00    109    M  =  391  ,0-30 
XL<M>=.0995 
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XWlMli.  0  24875' 

109  CONTINUE 

DO  110  M=UJ1  ,0.74 
XLCM)=.02*87S 
XW(M)=.024ft7S 

110  CONTINUE 

DO  111  M  =  475.4a2 
XUfM>=. 024875 
XWfM)=. 024975 

111  CONTINUE 

00    10    1=1, NP 

00    15    J=U10 

IF(M. NE.NE)   THEN 

IF(J.Eq.lO)    THEN 

REA0(5,«S)    SIGY,SIC7_,S1GX,SIGYZ 
ELSE 

R£AD<5,35)    3IGY,SIGZ,SIGx,oIGY2 
END    IF 
ELSE 

READ(5,30,I     SIGY,SIGZ,SIGX.SIGYZ 
GO   TO   SS 
£N0    IF 

w«Ur£(6,<+0)    M,XLfM)»XW(M>,SI6f,SIGZ,3IGX,SIGY2 
t5  CONTINUE 

10    CONTINUE 

53    *RJT£C6,40>    M.XL(W)  ,XW(H)  ,SrGY,SIGZ,SIGX,5IGrZ. 
30    F0ftMArci8X,O£l 1.0) 
35    F0RMAT(l3X,4E1l .0,//'/) 
'10    F0RMAT(3X,  I5,2F11  .7.4E12.4) 
45    FORMAT  (  18X,4£<1  .0  ,////////////  ) 

50    FQRMATflX, 'ELEMENT  NO . * ,4X , 'L* , J  OX,  '*• , 7X , • SIGY ' ,3X , • SIGZ ' 
*,8X,  'S1GX'  ,7%,  'SIGYZ'  ./) 
STOP 
ENO 
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